jyrjr 



biO 



PY 

noi 

»_ 

_\_ ~ 



- 0*1 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 





N / 5‘ I'l' 



AERDTHERMODYNAMICS OF A JET ENGINE TEST 
CELL FACILITY 

by 

Eric A. Nicolaus 

• * i* 

Septanber 1988 



Thesis Advisor: David Salinas 



Approved for public release; distribution is unlimited. 



T242206 



ECURITY 1 CLASSIFICATION OF THIS PAGE 



REPORT DOCUMENTATION PAGE 



a REPORT SECURITY CLASSIFICATION 

UNCLASSIFIED 



ID RESTRICTIVE MARKINGS 
NONE 



SECURITY CLASSIFICATION AUTHORITY 



3 DISTRIBUTION /AVAILABILITY OF REPORT 



D DECLASSIFICATION / DOWNGRADING SCHEDULE 



Approved for public release; distribution 
is unlimited. 



PERFORMING ORGANIZATION REPORT NUMBER(S) 



$ MONITORING ORGANIZATION REPORT NUMBER(S) 



a NAME OF PERFORMING ORGANIZATION 

Naval Postgraduate School 



60 OFFICE SYMBOL 
(If applicable) 

69 



7a NAME OF MONITORING ORGANIZATION 

Naval Postgraduate School 



c. ADDRESS (Cry, Srare, and ZIP Code) 

Monterey, California 93943-5000 



7d ADDRESS (City. State, and ZIP Code) 

Monterey, California 93943-5000 



a. NAME OF FUNDING / SPONSORING 
ORGANIZATION 



80 OFFICE SYMBOL 
(If applicable ) 



9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 



ADDRESS (City, State, and ZIP Code) 



10 SOURCE OF FUNDING NUMBERS 



PROGRAM 


PROJECT 


TASK 


WORK UNIT 


ELEMENT NO 


NO 


NO 


ACCESSION NO 



TITLE (Include Security Classification) 

AEROLHERMODYNAMICS OF A JET ENGINE TEST CELL FACILITY 



PERSONAL AUTHOR(S) 



NICOLAUS, ERIC A. 



la. TYPE OF REPORT 


13b TiME COVERED 


14 DATE OF REPORT (Year, Month, Day) 


15 PAGE COUNT 


Master's Thesis 


FROM TO 


1988 September 


151 



SUPPLEMENTARY NOTATION 



The views expressed in this thesis are those of the author and do not 
reflect the official policy or position of the Department of Defense or the U.S. Government. 





COSAT l 


CODES 


IB SUBJECT TERMS ( Continue on reverse if necessary and identify by block number) 


FIELD 


GROUP 


SUB-GROUP 


PHOENICS , Hush House, Test Cell, Aerothermal properties 



) ABSTRACT ( Continue on reverse if necessary and identify by block number) 



his thesis consists of a three-dimensional numerical analysis of the Jet 
ircraft Hush House located at Naval Air Station Jacksonville, Florida, 
tilizing the PHOENICS Code allows for the determination of the aerothermal 
haracteristics including velocity, pressure, enthalpy, turbulent kinetic 
nergy and the dissipation rate of turbulent kinetic energy in the facility 
uring testing of the U.S. Navy's F-4 (Phantom II) J-79-GE-8 gas turbine 
ngine with afterburner. How and by what method PHOENICS arrives at this^ 
olution is discussed. Of greatest importance is the resulting "behavior" 
f the aerothermal system. 

roblems encountered using the PHOENICS Code, resulting numerical solutions 
o the particular facility, comparison to actual test data and recommendations 
or further applications of the PHOENICS Code are presented. 



j distribution/ availability of abstract 
□ UNCLASSIFIED/UNLIMITED □ SAME AS RPT 


□ DTiC USERS 


21 ABSTRACT SECURITY CLASSIFICATION 

UNCLASSIFIED 


a NAME OF RESPONSIBLE INDIVIDUAL 

DAVID SALINAS 


22b TELEPHONE (Include Area Code) 

408-646-3426 


22c OFFICE SYMBOL 

Code 69Sa 



ID FORM 1473, 84 MAR 83 APR eot.on may De useo until exnausted SECURITY CLASSIFICATION OF THIS PAGE 



All otner editions are obsolete 



ft US Government Pfintinf Office It!*— *04 24 3 



1 



UNCLASSIFIED 



Approved for public release; distribution is unlimited 



Aerothermodynamics of a Jet 
Engine Test Cell 



by 



Eric A. -Nicolaus 

Lieutenant Commander, United States Coast Guard 
B.S., United States Merchant Marine Academy, 1978 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN MECHANICAL ENGINEERING 

from the 



NAVAL POSTGRADUATE SCHOOL 

September 1988 






ABSTRACT 



This thesis consists of a three-dimensional numerical 
analysis of the Jet Aircraft Hush House located at Naval Air 
Station Jacksonville, Florida. Utilizing the PHOENICS Code 
allows for the determination of the aerothermal 
characteristics including velocity, pressure, enthalpy, 
turbulent kinetic energy and the dissipation rate of turbulent 
kinetic energy in the facility during testing of the U.S. 
Navy's F-4 (Phantom II) J-79-GE-8 gas turbine engine with 
afterburner. How and by what method PHOENICS arrives at this 
solution is discussed. Of greatest importance is the resulting 
"behavior" of the aerothermal system. 

Problems encountered using the PHOENICS Code, resulting 
numerical solutions to the particular facility, comparison to 
actual test data and recommendations for further applications 
of the PHOENICS Code are presented. 
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New and overhauled aircraft gas turbine engines undergo 
stringent performance testing before returning to operational 
use. To allow for safe and effective monitoring of the engine 
and aircraft, gas turbine Test Cell and Hush House facilities 
are utilized to perform monitoring of key temperatures, 
turbine speeds and flow rates associated with the engine. The 
test facility additionally provides distortion-free air flow 
to the engine inlet (eliminating rain and wind effects), for 
lower exhaust temperatures exiting from the facility and 
reduced noise pollution to the surrounding area. Furthermore, 
these facilities allow for final engine "trim" settings to 
ensure that maximum power is available. 

Currently the test facilities designed and built for the 
Navy vary widely. The various facilities include: 

1. The Run-up pad--where a tied down aircraft conducts 
power checks in an open air environment. 

2. Sound suppressor--located at the Run-up pad to reduce 
noise . 

3. Open Test Stand (OTS) — located outdoors that supports an 
aircraft engine and has portable instrumentation 
attached. 

4. Hush House — a totally enclosed facility allowing around 
the clock testing of aircraft. 

5. Test Cell (TC)--an enclosed facility allowing testing of 
various engine sizes out of the airframe. 
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Design engineers have used experience and common sense to 
build these facilities with varying degrees of success. 
Aerothermal monitoring supplements the design of the 
structures . 

With the tremendous technological advances made in recent 
years, much higher temperatures and thrusts are available 
compared to similar size and weight engines of just a few 
years ago. Considering that nearly 67% of the Navy Test Cells 
and Hush Houses are at least 25 years old it is evident that 
a new generation of standardized facilities meeting federal, 
state and local noise, pollution, security and safety 
standards is required. [Ref.l:p.3] 

The Navy Environmental Health Center cited total claims 
against the government in 1981 at $114 million due to physical 



injury, 


hearing, or 


respiratory 


damage to 


military 


and 


civilian 


personnel . 


Additionally , 


the 1981 


ruling by 


the 



federal court system that military gas turbine test facilities 
must comply with Section 233 of the Clean Air Act accelerated 
the need for modern test facilities meeting the above 
requirements plus those projected for the future. [Ref.l:p.3] 
The goal of the Navy is to provide for standardized 
aircraft and engine test facilities that meet both safety and 
environmental concerns while providing the latest technology 
to monitor the engines. To accomplish this, complete 
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development of test facilities with reliable operation and 
maintenance throughout the facility life cycle is necessary. 
[Ref . 1 :p. 1] 

When contemplating the inherent design criteria for the 
new facilities, problems such as: turbulence, non-concentric 
flow patterns, multiple jet engine aircraft (including varying 
temperature and thrust arrangements), supersonic velocities, 
extreme temperature and pressure gradients, and awkward 
geometries must be appraised. [Ref. 5] 

Enigmas associated with the design of the structure may 
not be obvious. Complete information dealing with the 
aerothermal properties is a must in basing decisions on new 
multi-million dollar expenditures. Therefore, to minimize the 
possibility of an erroneous design, the PHOENICS Code is 
employed first on documented Test Cell and Hush House 
structures having actual field test data. After verifying the 
accuracy of PHOENICS with the test data, the Code eventually 
could be used at a Navy installation (such as the Naval Civil 
Engineering Laboratory) to design future test facilities 
providing accurate modeling data. 

In this instance, the performance of a Navy J-79-GE-8 gas 
turbine engine located at Naval Air Station Jacksonville, 
Florida Hush House facility is investigated. The data is used 
to determine the behavior of the existing structure and will 
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be employed to assist engineers in the research 
development (R&D) phase of the facility design project. 



and 
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II . PROBLEM FORMULATION 



A- HUSH HOUSE CHARACTERIZATION 

The Hush House under investigation is an "Air Force" type 
located at Naval Air Station Jacksonville, Florida. In 1980 
the Air Force purchased 25 and later 54 demountable Hush 
Houses with Test Cell capability. This standardization 
accompanied with such a large scale purchase and the 
flexibility of aircraft or engine test capabilities allowed 
for a cost effective combined design. Because of the large 
scale purchase the cost per Hush House was a relatively low 
$2.2 million. The only major modification performed by the 
Navy to the Air Force design was inclusion of an additional 
AFFF fire fighting system. The Hush House design does not 
address air pollution standards. Some type of air pollution 
abatement system will undoubtably be required in newer 
designs. [Ref.l:p6] 

The general dimensions of the Hush House is 28.7 meters 
(94 feet) wide by 25.5 meters (83.7 feet) long with the oval 
shaped "augmenter tube" extending an additional 31.3 meters 
(103 feet) from the rear of the Hush House (See Figures 2.1 
and 2.2). The aircraft is placed so that the engine exhaust 
is approximately 4.6 meters (15 feet) from the augmenter tube 
entrance. The most common construction material for the 
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structure is steel. The design is geometrically simple 
allowing for mass production. The air flow in the last two 
inlets nearest the augmenter tube provides both engine intake 
and cooling air for the surrounding engine area. The final 7 
meters (23 feet) of augmenter tube is angled at approximately 
45 degrees running from the floor upward to provide for a 
smoother directional transition of gases. Acoustic damping 
materials line the augmenter tube, inlets and aircraft bay 
walls to reduce exhaust noise. 

The Hush House is used in support of the Naval Air Rework 
Facility (NARF) located at the Naval Air Station. Although the 
Hush House is presently utilized for engines fully mounted in 
an airframe, it is anticipated that engine "Test Cell" testing 
(where an engine is mounted on a test stand) will be approved 
in the future. The facility is normally in operation from 0700 
to 2400 Monday through Friday with an average of 1 to 5 
aircraft being tested per day. One of the most common aircraft 
tested is the F/A-18 Hornet powered by two General Electric 
F-404 gas turbine engines. Typically, after completion of 
overhaul of the aircraft and engines, the aircraft will be 
towed to the facility and backed into the Hush House. The hugh 
aircraft inlet doors to the facility will then be closed to 
provide an enclosed area free of foreign matter, outside wind, 
rain, and other environmental and personnel concerns. 
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With the aircraft fully enclosed and securely tied down, 
a normal crew of three mechanics will perform final engine 
"trim". One person is located in the Hush House Control Booth 
to ensure manning of the fire fighting systems. The second 
will physically adjust the fuel control mechanism on the 
engine via radio signals received from the cockpit where the 
third crewman is. The engines are controlled during testing 
by the third crewman in the aircraft who is additionally 
responsible for monitoring the various fuel, oil, % power, 
hydraulics, and other gauges. Testing normally lasts for a 
period of two hours. Only one engine is tested at a time with 
the other normally at idle. Air intake and cooling is provided 
by five inlet areas located on either side of the Hush House 
(See Figures 3.2 and 3.3). Wire screening, air diffuser and 
noise abatement equipment is present at the inlets. Installed 
fire fighting and wash-down equipment is present in the 
aircraft bay area with the air inlets capable of being 
automatically closed. 

B. ISSUES CONCERNING THE HUSH HOUSE 

The Hush House is currently performing as designed with 
the below exceptions. Investigation indicates some cracking 
of the welds which fuse together the stainless steel 
"corrugated" sheet metal overlaying the sound absorption 
material located in the augmenter tube. Also, actuation of 
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heat sensors in the overhead of the aircraft bay sometimes 
occurs when the engines are idled for a period of time. 
Apparently some of the exhaust gases travel upward when low 
mass flow rates are present instead of going down the 31 plus 
meters (103 feet) of augmenter tube. 

In general, personnel at the Naval Air Station were 
pleased with the overall operation of the Hush House. Although 
vibration of the entire physical structure has been noted in 
this facility as in other designs, it does not appear to cause 
any structural problems such as the requirement for inspection 
and retightening of bolts as occurs in other Navy Hush House 
designs . 

Physical test data is available for an identical Air Force 
facility for the F-4 Phantom II and is included in parts of 
this report for comparison with results from the PHOENICS 
Code. The test data provides an important comparison for this 
thesis to verify the computer model. [Ref. 2] 

A physical inspection of an identical Hush House to Naval 
Air Station Jacksonville's was under-taken in June 1988 at 
Dobbins Air Force Base in Atlanta Georgia. Discussions with 
the various Air Force personnel corroborated the previously 
mentioned problems with the exception of the heat sensor 
activation problem. The Dobbins Air Force Base Hush House was 
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utilized for both Test Cell and Hush House operations 
involving the F-15 (EAGLE) fighter. 
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Air Inlets 




Figure 2.1 Actual Hush House (Top View Looking Down) 
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Aircraft Entrance Doors 



Figure 2.2 Actual Hush House (Side View) 
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XXI. THE PHOENICS CODE 



A - AIM EXPLANATION 

PHOENICS, which stands for "Parabolic Hyperbolic Or 
Elliptic Numerical Integration Code Series", is a computer 
code capable of simulating fluid flow (including changes of 
state), heat transfer and chemical reactions. The code is 
copyrighted by CHAM (Concentration Heat and Momentum) 
Limited which was founded by D. Brian Spalding in 1972. 
[Ref . 3:p. 1 . 1] 

As shown in Figure 3.1, the PHOENICS Code is structured 
into three major components consisting of "EARTH", 
"SATELLITE", and "GROUND STATION". The "EARTH" program may be 
considered as a "black box" which contains the computational 
procedures for all applications. The "SATELLITE" program 
consists of several data-input methods which define the 
boundaries and scope of the problem. Finally, the "GROUND 
STATION" program has various subroutines containing built- 
in physical models which may or may not be selected depending 
on the program application. 

By default "EARTH" can solve for up to 25 dependent 
variables. In the Hush House application these include 
pressure, all three dimensional velocity components, 
enthalpy (temperature), turbulent kinetic energy, and the 
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turbulent kinetic energy dissipation rate. The dependent 
variables can be enlarged if desired. 




Figure 3.1 Satellite, Earth and Ground 

The PHOENICS Code is exceptionally flexible allowing for 
varying sizes and scales of structures from "micro" to 
"macro". A "final" solution is highly dependent on the 
accuracy required by the user as well as the grid sizing used. 
The higher the accuracy and/or number of "cells" desired, the 
greater the computational effort required for a solution. In 
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this investigation, a "grid independent solution", where a 
sequence of "finer" grids would be used, was not possible due 
to time limitations. 

In order to solve a given problem the user must create a 
"Ql" file. The "Ql" file describes all aspects of the system 
being investigated. This file is partitioned into several 
groups allowing the user to specify various aspects of the 
physical system. Examples are structure size and shape. 
Additionally, other groups are capable of providing the 
necessary grid spacing, boundary conditions, engine mass 
flow rate, engine exhaust temperature and printing options. 
The "Ql" file is read by the "SATELLITE" program which is 
recognized and processed by "EARTH". The "EARTH" program 
solves discretized conservation equations. A typical 
conservation equation (before discretization) which "EARTH" 
solves can be expressed in the form: 

6 

-- (rp<f>) + div(rpV$ - rT* grad $) = rS* 

5t 

where : 

r = phase volume fraction 

<J> = dependent variable 

T* s exchange coefficient (laminar or turbulent) 

S» 2 source or sink term 
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p s density 
V = velocity vector 

Here $ is the dependent variable currently in question. For 
example, using the conservation of energy equation, <t> would 
be enthalpy and T* is the diffusion constant. [Ref . 3 : p . 2 . 4 ] 
The values computed and displayed for each cell are the 
solutions of the algebraic equations relevant to the laws of 
physics for that cell. The laws of physics are those of 
conservation of mass , energy and momentum and any other 
dependent variable. The general form of the algebraic 
linearized equations resulting from a discretization of the 
previous field equation is: 

aE $e + aw<J>w + aN<J>r. + as$s + aH<J>H + aL'J’L + ai$i + S 
= 

aE + aw + as + as + an + ai + ai + ap 

where $ f > is the value of the dependent variable at a 
particular cell and subscripts E, W, N, S, H and L are 
neighboring locations. The values of <t> at these locations are 
known from a previous iteration. Quantities S and a P express 
the influence of a source of the entity <t>. The a's are 
influence coefficients, temporarily treated as constants, with 
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the above subscripts expressing interaction between neighbor 
cells. [Ref .3:p.2.4] 

Solving the multitude of algebraic equations is 
complicated and time consuming. An iterative manner is used 
until the imbalance between the left and right sides of all 
the above equations are negligible. Several schemes can be 
used to perform the "sweeps" (where a sweep solves for the 
dependent variables for each cell in the grid one time) 
necessary for the iteration process. The results will be 
identical although the time to complete the task will vary 
with the iteration scheme. [Ref . 3 : p . 2 . 5] 

There are several ways to determine whether the problem 
has "converged" to a solution. Results are monitored to 
determine whether values of dependent variables are remaining 
constant or fluctuating. The "error" (residuals) is tabulated 
for each sweep at the end of the computer printout. By 
observing "residuals" from each sweep a pattern of decreasing 
("converging") values should be present. If the "residuals" 
are increasing ("diverging") adjustments of several constants 
which in-turn apply a dampening effect upon the variables in 
question may be required. The printout of successive residual 
values demonstrate convergence by a reduction in the 
residuals. That is, a visually downward sloping line on the 
accompanying graph of residuals shows convergence. 
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Conservation balances on each dependent variable can be 
verified using "net source" values which are tabulated at the 
end of each run. The mass flow rate (mass flow rate = density 
* area * velocity) is additionally calculated, in this 
instance, for the five inlets and one outlet confirming that 
no hidden, missing, or improper boundaries are present. The 
problem has "converged" when an acceptable level of error 
remains . 

Successful "convergence" of the PHOENICS program is 
dependent upon experience and an understanding of the results. 
If the residuals of a dependent variable show divergence, the 
user can manipulate the rate of change of values of that 
variable by using the RELAX command to return the variable 
back to the convergence path. Relaxation is placed on a 
variable to control the variables maximum change from each 
sweep. This minimizes the possibility of "over-correction" 
or "blowing up" of the values. Smaller values of relaxation 
result in smaller value differences per iteration. However, 
excessively small relaxation values require additional 
computer time. A balance gained by user experience enables 
optimal convergence with minimal computer time. 

When multiple runs are required for large programs output 
data is stored in a "DF09" file which allows the next run to 
resume the computations (i.e. f a restart) where they last 
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ended. The complexity of the problem and the required 
accuracy determine the quantity of computer runs needed. The 
number of runs can vary from the tens to the thousands and 
higher . 

B . ASSUMPTIONS IN THE MODEL 

There are several approximations that deserve mention in 
the modeling of the Naval Air Station, Jacksonville, Florida 
Hush House. Because the geometry of the model uses cartesian 
coordinates, the actual Hush House's rounded and curved 
surfaces, including the augmenter tube, Hush House overhead, 
engine and exhaust gas exit ramp, were approximated as flat 
stepped surfaces. 

The inlet "muffler" arrangement has been eliminated for 
simplification. The corrugated metal found throughout the 
structure to cover the sound absorption material was 
considered to have a friction surface. The model, therefore, 
included friction for all walls, floors, etc. 

The actual test data used for comparison between the 
model and facility contain several notable differences. The 
model uses the symmetry of the centerline of the Hush House 
to reduce the number of cells by half. The engine is located 
on the centerline in an ideal cradle position. Because of 
the symmetry only half of the engine is observed as if the 
Hush House was being used as a test cell without the aircraft 
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present. Since the correct engine exhaust area is used the 
units of mass flow rate and temperature are correct for the 
model configuration. The engine is configured for afterburner 
operation. The physical aircraft is not present. 

Actual test data used for comparison consists of the 
entire aircraft being placed in the facility. The engine 
configuration has both engines running with one at 
afterburner while the other engine is at "engine idle". The 
obvious result of the engine arrangement is a non-symmetr ical 
engine location in comparison to the symmetric model. 

C. MODELING THE PROBLEM 

Data, including blueprints of the Hush House and aircraft 
engine information, was gathered to accurately depict the 
facility configuration. The information obtained was in 
English units and was converted to the International Metric 
System or "S.I." using standard conversion factors. 

With the appropriate information gathered it was time to 
"fit" the Hush House to the cartesian coordinate system of X, 
Y and Z coordinates using PHOENICS version 1.3. The cartesian 
coordinate system was choosen because of its simplicity. Since 
the actual structure does have rounded surfaces some 
approximations were made specifically to the augmenter tube 
and overhead of the inlet bay. Additionally, the aircraft was 
eliminated to reduce the complexity of the problem allowing 
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the structure to serve in its "Test Cell" configuration. 
Drawings are found in Figures 3.2 and 3.3 of this section 
illustrating the key items found in the structure. The 
drawings view through the structure to note the key items, 
such as: air inlets, augmenter tube and engine locations. 
Additionally, the structure walls were considered to be of 
infinitesimal thickness to further simplify the grid 
structure. The engine was simplified from a cylindrical shape 
to one of a rectangle of equal intake and exhaust area as the 
original cylinder. Other small items, such as the engine 
stand, were eliminated. 

To take advantage of the symmetry of the Hush House and 
thereby reduce computer CPU time, only half of the Hush House, 
engine and augmenter tube were modelled in the "Ql" file. This 
was easily accomplished by dividing the structure at the 
centerline. This also appreciably reduces the size of the 
"Ql" file. 

Figures 3.2 and 3.3 graphically display the X-Z and the 
Y-Z planes of the structure. Intuition, tr ial-and-error , and 
available computer time dictated the grid spacing used. 
Obviously, volume areas of interest such as the engine exhaust 
and augmenter required smaller cell size. 

The "Ql" file describes the structure using several types 
of statements. The "CONPOR" statement allows the user to 
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define walls, overheads and the floor. The "PATCH” and 
"COVAL" statements work together to indicate items such as 
an engine exhaust or atmospheric boundaries. Dependent 
variable values for say the atmospheric boundaries can be 
"fixed" or allowed to "float" depending on the problem 
definition. When the walls, inlet bay overhead, augmenter 
tube, engine and floor dimensions have been measured, a grid 
for the system can be constructed, based on locations and 
anticipated responses over different "areas" of the system. 
Areas of anticipated large response gradients would have 
smaller cell volumes than low activity areas. 

Note that in figure 3.2 the Z-axis consists of 40 cells 
consecutively numbered from 1 to 40 going left to right. The 
X-axis consists of 10 cells consecutively numbered 1 through 
10 going from lower to upper. Figure 3.3 is similarly numbered 
where the Y-axis runs from; 1 to 12 going from lower tc upper 
and the Z-axis is the same as previously described. 

In addition to describing the geometric structure and 
prescribing the "P.ELAXATION" devices, the properties of the 
fluid medium must be specified. The following boundary 
conditions were used: 

1. ambient temperature = 295 K 

2. ambient pressure = 101325 Pa 

3. gas constant = 286 Joules/(kg - K) 
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4. Cp = 1004 Watt-sec/(kg - K) 

5. mass flow rate (engine) = 197.72 kg/(sec-sq meters) 

6. temperature (engine out) = 1944 K 

Additionally, "friction" was applied to all surfaces of 
the Hush House. Appendix A contains the final form of the 
Q1 file for reference. An in-depth explanation of all the "Ql" 
program parameters can be found in the "Beginners Guide and 
Reference Manual" (See Reference 4). 

Several methods and theories have been previously used to 
attempt to reach a "converged" solution in other test 
facilities. The method chosen for this project was to 
initially minimize the "excitation" of the model thereby 
reducing the "jump" required in solving the momentum, 
enthalpy, turbulent kinetic energy dissipation rate, turbulent 
kinetic energy and mass equations. To accomplish this, during 
successive runs (where a "run" consists of several "sweeps"), 
small incremental increases to the source were applied, i . e 7 
jet engine exhaust temperature and mass flow rate were slowly 
increased, from, run to run, while carefully monitoring the 
residuals as well as the mass flow rate in and out of the 
structure. The turbulence model, "KEMODL" , was used 
throughout rather than considering the flow first as laminar 
and then adding turbulence. This method appears to have met 
with considerable success. The table that follows illustrates 
the iterative process used throughout the computer modeling: 
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This amounts tc approximately 147 hours of computer CPU time 
on the IBM 3033 mainframe computer installed at NPS. 

Using the MVS "batch" system on the IBM 3033 mainframe 
proved to be time consuming. The PHOENICS Code is designed to 
work both interactively and in batch mode. Due to the high 
computer usage at NPS the "interactive" mode was not 
available. Limited to the batch mode, careful planning for 
evening or night runs was required to maximize CPU time. The 
number of "Degrees of Freedom" (DOF) that must be solved for 
gives an idea of the problem complexity. DOF = NX * NY * NZ 
* the 7 dependent variables. Solving this for the grid size 
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of 10 * 12 * 41 yields nearly 35,000 values calculated per 
each sweep. Typically, 150 to 200 sweeps could be accomplished 
during the day using four 15 minute (class "G") runs of 50 
sweeps each. After normal working hours three one hour (class 
”j") runs totalling 500 to 600 sweeps were accomplished. Data 
stored on the all important DF09 file was continuously 
"backed-up" after a successful run to prevent loss of data if 
an improper subsequent run resulted in erroneous data. The 
backed-up data could then be restored to the DF09 file to 
enable continuation of the runs without starting at the very 
beginning . 
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Figure 3.2 Hush House Model (X-Z plane) 
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Figure 3.3 Hush House Model (Y — Z plane ) 
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XV - RESULTS 



A- THE FINAL SOLUTION 

After obtaining the "solution" for the problem the 
question still remaining is just how "good" are the computer 
model results? Some of the "tools" available for use, such as 
checks on residual values and mass flow rate calculations, 
were discussed in Chapter III. The checks indicate that a 
converged solution is obtained, but, still do not prove the 
computer model data validity with the actual facility data. 

Using actual (experimental field measurements) results 
from Reference 2 allows for an accurate evaluation of the 
computer model. Appendices "C" and "D" offer some experimental 
data for velocities and temperatures, thereby indicating mass 
flow rate and augmentation ratio. Additionally, figures 4.1, 
4.2 and 4.3 found in this section give the numerical 
(computer) model results of differential pressure, temperature 
and axial velocity profiles encountered at the IX=1 slab. 
Graphics generated pictures (generated by "PHOTON" which 
converts the computer generated results to graphical format) 
allow for the interpretation and digestion of the nearly 300 
pages of numerical output of a PHOENICS run. A condensed 
output is found in Appendix "B". In this case, PHOTON graphics 
specifically provides graphical information for a better 
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interpretation of velocity fields and both pressure and 
temperature profiles and contours enabling a better 
understanding of computer generated data. Use of all these 
additional "tools" allow for a determination as to whether the 
computer results are valid. 

Observing figures 4.4 and 4.5 (top looking down) show the 
flow pattern of the air entering the five air inlets into the 
Hush House at IY=4 (engine level). A "typical" low velocity 
air flow region (Figure 4.4), illustrated by the minuscule 
length of the velocity arrows (velocities of less than 2 
meters/sec or 6.5 feet/sec), is noted between the first two 
and final three inlets. Figures 4.6, 4.7 and 4.8 illustrate 
the side view velocity flow patterns for the overall facility, 
engine area, and exhaust to atmosphere respectively. Note the 
ever reducing velocity magnitudes which vary from 818 meters 
per second (2684 feet/sec) at the engine exhaust to less than 
100 meters per second (328 feet/sec) as the gases travel 
further away from the engine and down the augmenter tube until 
reaching the exhaust. 

Figures 4.4, 4.5 and 4.9 (top looking down) show the air 
and gas flow patterns entering and exiting the engine region 
and augmenter tube at engine level. Figure 4.10 indicates the 
flow patterns at the top of the augmenter tube (IY=8) for the 
entire Hush House. Of interest is the counter-flowing currents 
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at the top of the augmenter tube starting at the engine 
exhaust and running to nearly half way down the augmenter tube 
itself . 

The engine occupies cells defined by IX=1, IY=4 to 5 and 
IZ=12 to 14. As anticipated, figures 4.4, 4.9 and 4.11 show 
that at engine level the air in the Hush House moves toward 
the engine inlet from the surrounding area while air also is 
entrained at the exhaust of the engine. Figure 4.1 illustrates 
this entrained "effect" noted by the high differential 
pressure of the exhaust gases immediately exiting the engine 
at IX=1 and IY=4 . Less than a meter after exiting the engine 
exhaust, the exhaust gases and inlet air have begun mixing 
which results in an area of drastic temperature reduction for 
the exhaust gases in the Hush House and immediate augmenter 
area. Gradually, as the gases pass further downstream, 
pressure begins rising which reduces the velocity of the gases 
in the augmenter tube until the gases exit to the atmosphere. 
At the inlet area to the engine (at Z=18 meters on figure 4.1) 
a low pressure area is noted as one would also anticipate as 
the air is drawn into the engine inlet. 

The exhaust gases exit the engine at 1944 Kelvin (3040 F) . 
The inlet air and exhaust gases mix due to the entrainment 
"effect" causing a distinct drop (See figure 4.2) in 
temperature as the gases travel further down the augmenter 
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tube. Temperature contours at the IX=1 slab shown in figures 
4.12 and 4.13 portray the temperature patterns out of the 
engine progressing down and out of the augmenter tube. 
Analogous temperature contours at the IX=3 slab are presented 
in figures 4.14 and 4.15. Figures 4.16, 4.17 and 4.18 embody 
the temperature contours at the IY=5 (engine) level while 
figure 4.19 shows the temperatures at the top of the augmenter 
tube at IY=8. The temperature contours indicate a 
progressively widening area of decreasing exhaust gas 
temperatures where the vast majority of the temperature drop 
occurs in the first five meters after the engine exhaust (the 
temperature differential goes from 1944 K to 800 K) . 

Differential pressure contours of the engine exhaust at 
IX=1 are shown in figures 4.20 (IZ= 14 to 24) and 4.21 (IZ= 
25 to 40) . As indicated in the pressure profile, the exhaust 
jumps to nearly 6000 pascals above ambient upon exiting the 
engine, but soon after drops into a low pressure region until 
reaching ambient pressure half way down the augmenter tube at 
IZ=24. Gradually, the pressure builds until the gases exhaust 
to the atmosphere. Similar contours for the IX=3 slab are 
included in figures 4.22 and 4.23. Pressure contours for the 
X-Z planes at IY=5 (engine level) and IY=8 (top of the 
augmenter tube) are contained in figures 4.24, 4.25 and 4.26. 
Of noted interest is how the differential pressures remain 
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fairly uniform at a given IZ location as the gases travel down 
the augmenter tube with the exception of the region initially 
downstream from the engine exhaust. 

B . ACTUAL VERSUS MODEL. COMPARISONS 

Generally, the results from the computer model closely 
approximated the actual. The below "actual" data, and that of 
Appendix "C" and "D" for the computer model are indicated in 
both "English" and "S.I." units below. A comparison of the 
PHOENICS output visually shown in Appendices "C" and "D" 
illustrates the dissimilarity between the actual and computer 
model. As indicated by the data, the differences in the 
computer model (235.5 ft/sec or 71.8 meters/sec) and actual 
(227.2 ft/sec or 69.3 meters/sec) averaged velocities found 
in the augmenter tube at slab IZ=27 differed by less than 4 
percent. Averaged temperatures in the same location differed 
between 401.6 degrees Fahrenheit (478 K) for the computer 
model and 454 degrees Fahrenheit (507 K) actual with a 
resultant difference of 6 percent. Reference 2 does not 
indicate whether the effects of radiant energy from the engine 
exhaust are considered for the actual temperature readings 
taken in the augmenter tube. If these effects were not 
considered the 6 percent difference would be reduced. The Hush 
House mass flow rate for the computer model (2373 lbm/sec or 
1076 kg/sec) versus the actual (2067 lbm/sec or 938 kg/sec) 
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differed by 14.8 percent. The computer model mass flow rate 
value was actually doubled (hence the 2373 lbm/sec or 1076 
kg/sec value) to "model" the entire engine and Hush House. 
The Augmentation Ratio is defined as follows: 

M Inlet — Me n g 1 n e 

A. R . = 

Mi nlet 



where : 

Mi n 1 e t = inlet mass flow rate 

Mengine = engine mass flow rate 

The inlet mass flow rate = air density * five inlet 
opening areas * averaged inlet velocity. The engine mass flow 
rate is given at 197.72 kg/sec/sq meters * 0.195 sq meters 

engine exhaust area for half of the engine. 

Augmentation ratios range from 5 to 15 depending on the 
design of the Hush House or Test Cell where the higher the 
ratio typically indicates a more efficient design. The 
computer model augmentation ratio was found to be 10.6 while 
the actual was 8.3. The resultant difference is therefore 
nearly 28 percent. Considering the limited data taken at this 
actual power level (See Appendix C) it is quite conceivable 
that the actual augmentation ratio could be greater. 

For the actual results the velocity, pressure and 
temperature measuring devices were photographed because of the 
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large fluctuations of the gauges. Additionally, velocity data 
at the final three inlets could not accurately be taken making 
it difficult to accurately determine the augmentation ratio 
and the true velocities at the inlets. 
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Figure 4.1 Differential Pressure Profile at Volume Slab 
IX=1 for IY=4 (Engine Height) 
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Figure 4.2 Temperature Profile at Volume Slab IX=1 
for IY=4 (Engine Height) 
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Figure 4.3 Axial Velocity Profile at Volume Slab IX=1 

for IY=4 (Engine Height) 
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Figure 4.4 Flow Vector Plot of Hush House Inlets and 
Engine Area (X-Z Plane at IY=4) 
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Figure 4.5 Flow Vector 
{Engine Exhaust 



Plot in X-Z Plane at 
and Augmenter Tube) 



IY=4 
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Figure 4.6 Flow Vector Plot in Y-Z Plane at IX=1 

(Entire Grid) 
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Figure 4,7 Flow Vector Plot in Y-Z Plane at IX=1 

(Engine Exhaust Area) 
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Figure 4.8 Flow Vector Plot in Y-Z Plane at IX=1 
(Exhaust to Atmosphere) 
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Figure 4.9 Flow Vector Plot in X-Z Plane at IY=4 
(End of Augmenter Tube) 
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Figure 4.10 Flow Vector Plot in X-Z Plane at IY=8 
(Mirrored Image of Entire Grid) 
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Figure 4.11 Flow Vector Plot of Augmenter Tube Area 

for X-Z Plane at IY=4 
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Figure 4.12 Temperature Contour in Y-Z Plane at IX-1 
(Engine and Augmenter Tube) 
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Figure 4.13 Temperature Contour in Y-Z Plane at IX=1 
(Augmenter Tube and Exit) 
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Figure 4.14 Temperature Contour in Y-Z Plane at IX=3 

(Engine Exhaust) 
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Figure 4.15 Temperature Contour in Y-Z Plane at IX=3 
(Engine Exhaust to Atmosphere) 
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Figure 4.16 Temperature Contour in X-Z Plane at IY=5 
(Mirrored Image of Engine Exhaust) 
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Figure 4.17 Temperature Contour in X-Z Plane at IY=5 

(Engine Exhaust) 
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Figure 4.18 Temperature Contour in X-Z Plane at IY=5 
(Engine Exhaust in Augmenter Tube) 
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Figure 4.19 Temperature Contour in X-Z Plane at IY=8 
(Engine Exhaust in Augmenter Tube) 
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Figure 4.20 Differential Pressure Contour in Y-Z Plane 

at IX=1 (Engine Exhaust) 
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Figure 4.21 Differential Pressure Contour in Y-Z Plane 
at IX=1 (Engine Exhaust to Atmosphere) 
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Figure 4.22 Differential Pressure Contour in Y-Z Plane 
at IX=3 (Engine Exhaust and Augmenter Tube) 
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Figure 4.23 Differential Pressure Contour in Y-Z Plane 
at IX=3 (Exhaust to Atmosphere) 
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Figure 4.24 Differential Pressure Contour in X-Z Plane 

at IY=5 (Engine Exhaust) 
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Figure 4.25 Differential Pressure Contour in X-Z Plane 
at IY=5 (Augmenter Tube) 



58 





Figure 4.26 



Differential Pressure Contour in X-Z Plane 
at IY=8 (Augmenter Tube) 
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V. SUMMARY AMD RECOMMENDATIONS 



A. SUMMARY 

Using the Parabolic, Hyperbolic, or Elliptic Numerical 
Integration Code Series (PHOENICS) a computer model of the 
Naval Air Station Jacksonville Florida's Hush House facility 
with a J-79-GE-8 gas turbine engine was simulated. The 
performance of this computer model was compared to actual 
performance data taken with an F-4 Phantom II undergoing 
testing at an identical facility located at Kelly AFB . An on- 
site inspection was conducted at Dobbins AFB in Atlanta 
Georgia to provide "hands-on" experience with the facility. 

The computer analysis consisted of slowly building up 
excitation in the model (using the KEMODL turbulence model) 
until full power was obtained. Computer generated (tabulated 
and graphical) results were compared with the actual test 
data. Results of the computer model, in comparison to the 
actual data, were in good agreement. Flow directions, 
velocities and temperatures were consistent with actual 
results . 

B . RECOMMENDATION S 

A comparison of actual to computer model data for this 
thesis indicates that PHOENICS is a useful design and 
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troubleshooting tool allowing the user to better understand 
the physics of the model. PHOENICS at NPS is now just coming 
into potent use with a variety of research projects for 
aeronautical and mechanical engineers. 

To ensure proper and full use of this code several major 
"glitches" must be overcome. Time constraints placed on the 
use of the IBM 3033 main-frame computer limit the number of 
"sweeps" obtainable by the student. The batch system does not 
lend itself to correcting possible errors that may arise. An 
interactive system, such as the "VAX" stand alone systems 
being installed in Mechanical Engineering, is an ideal answer. 

PHOTON was not available with the PHOENICS installation 
on the IBM 3033 main-frame computer. A complicated conversion 
process (See Appendix "E") was undertaken to convert IBM to 
VAX "format" to allow the use of PHOTON using data from the 
main-frame. This conversion process which involved a swap of 
information from IBM to DEC to VAX hardware took nearly 2 
weeks to accomplish. In September 1988 this transfer procedure 
will no longer be possible at NPS with the elimination of 
certain hardware at the Computer Science Department although 
data transfer will still be possible to other facilities such 
as CHAM. 

The use of PHOTON early in the computer "runs" would have 
uncovered several missing "CONPOR" statements and therefore 
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"verified" the structural integrity of the facility. 
Additionally, the amount of time required to "wade" through 
the 300 pages of printout could be reduced to several video 
screens with a much reduced volume of printout. 

The computer model used in this thesis was formulated 
using cartesian coordinates. Several parts of the structure, 
including the augmenter tube and engine, were curve shaped. 
Geometric approximations were made to substitute for the curve 
shaped areas. A more computationally time consuming model with 
increased numbers of cells and/or the use of body-fitted- 
coordinates (BFC) would enhance the computer model. 

As noted in Figure 4.2, over fifty percent of the 
reduction in engine exhaust gas temperatures occurred in the 
first 20 percent of the augmenter tube. Could an improved 
design include shortening the augmenter tube? This type of 
"what if" question can be answered using PHOENICS. With the 
variable data file ("DF09") already in place the grid pattern 
can be modified to shorten the augmenter tube. Minimal 
computations would be required to gather definitive results. 

Parameter studies, including shortening the length of the 
augmenter tube (to determine the effect of augmenter tube 
length on performance), movement of the engine location in the 
Hush House, changing the inlet locations and sizes, and the 
use of different turbulence models could lead to improvements 
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on the present Hush House design. In fact, the ability to 
change parameters and observe their effect on system 
performance is the prime attribute of numerical modeling. 

Finally, a discussion of turbulence is in order. PHOENICS 
offers four turbulence models which can be readily selected. 
The four models are: Algebraic, Mixing Length, K-Mixing Length 
or K-Epsilon. The user can also formulate a model of his own 
choosing which can be incorporated into the code via the 
"GROUND" program. Two questions arise. First, which turbulence 
model is appropriate for the problem under investigation? 
Secondly, how strong is the coupling between other system 
response variables and turbulence variables? The second 
question might be answered by running the PHOENICS Code on the 
problem using several different turbulence models. 
Unfortunately there was not adequate time to perform this in 
the investigation. The first question can not be resolved 
until a better understanding of turbulence is obtained. 

C . CONCLUS IONS 

The computer model results generated by this thesis 
closely parallel the actual test data. The results therefore 
indicate that the PHOENICS code can be successfully utilized 
to model future Test Cell and Hush House applications as well 
as "trouble-shoot" suspect areas of interest. Also, areas that 
may be difficult to properly monitor with test equipment could 
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be modelled instead. With the use of additional "cells" and/or 
BFC a more refined model could more closely duplicate the test 
facility providing more accurate results. Follow-up parametric 
studies to improve and verify the design and data seem in 
order. The results indicate that the Hush House is performing 
as intended. 
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APPEND XX A 
THE PHOENICS Q1 PILE 



TALK=F,RUN(1,1) 

* NAVAL AIR STATION JACKSONVILLE FLORIDA HUSH HOUSE TEST * 

* FACILITY. THESIS PROJECT OF LT ERIC A. NICOLAUS USCG * 

* WITH THESIS INSTRUCTOR PROF DAVID SALINAS. * 



* LIST OF VARIABLES * 

* * 

* WJET: VELOCITY OF ENGINE EXHAUST (METERS/SEC) * 

* MJET : MASS FLOW RATE OF GAS ( KG/SECOND/SQ METER) * 

* TJET : JET EXHAUST GAS TEMPERATURE (DEGREES KELVIN) * 

* P JET : EXHAUST PRESSURE OF JET ENGINE (NEHTONS/SQ METER)* 

* RHOJET: DENSITY OF THE GAS (KG/CUBIC METER) * 

* HJET i ENTHALPY (JOULES/KG OR NEHTON*METERS/KG ) * 

* KE JET : KINETIC ENERGY (JOULES OR NEWTON-METERS ) « 

* GASCON: GAS CONSTANT ( JOULES/l KG-DEG KELVIN)) * 

* TAMB: AMBIENT AIR TEMPERATURE (DEGREES KELVIN) * 

* HAMB: AMBIENT ENTHALPY (JOULES/KG) * 

* KEINJ : TURBULENT KINETIC ENERGY (JOULES) OF JET * 

* CSUBP : SP HT CONST PRE ( JOULES/( KG-DEG KELVIN)) * 

* EPINJ: DISSIPATION RATE OF TURBULENT KE (JOULES) OF JET* 

* KEINA: TURBULENT KINETIC ENERGY (JOULES) AMBIENT * 

* EPINA: DISSIPATION RATE OF TURBULENT KE (JOULES) AMB * 

* PI: PRESSURE (NEWTONS/SQUARE METER) * 

* HI: ENTHALPY (JOULES/KG) * 

* Ul: VELOCITY IN X-DIRECTION (METERS/SEC) * 

* VI: VELOCITY IN Y-DIRECTION (METERS/SEC) * 

* HI: VELOCITY IN 2-DIRECTION (METERS/SEC) * 

« TMP1 : TEMP OF THE FIRST PHASE (DEGREES KELVIN) * 

* RHOl : DENSITY OF THE FIRST PHASE (KG/CUBIC METER) * 

* PRESSO: REFERENCE ATMOSPHERIC PRESSURE ( NEWTONS/SQ mtr)* 

* TMP1A : GROUND REF TEMP USED IN GRND2 (DEG KELVIN) * 

* TMPIB : 1/CSUBP USED IN GRND2 EQU (KG-DEG KELVIN/JOULES)* 

* DRH1DP: DELTA RHO/DELTA PRESSURE H/RESPECT TO HI * 

* FACT: CORRECTION OR “FUDGE" FACTOR * 

* ENUL : LAMINAR KINEMATIC VISCOSITY (METERS SQ/SEC ) * 

* ENUT: TURBULENT KINEMATIC VISCOSITY (METERS SQ/SEC) * 

GROUP 1. RUN TITLE AND OTHER PRELIMINARIES 
TEXT! NAS JACKSONVILLE HUSH CHAMBER) 

REAL( WJET, MJET,TJET,PJET, RHOJET, HJET, KE JET, GASCON, CSUBP, TAMB) 
REAL( HAMB, KEINJ, EPINJ, KEINA, EPINA) 

GROUP 2. TRANSIENCE i TIME -STEP SPECIFICATION ( DEFAULT ^STEADY ) 

GROUP 3. X-DIRECTION GRID SPECIFICATION 
NX=10 

XF RAC ( 1 ) = -l j XFRAC ( 2 ) = 0 . 31 
XFRAC ( 3 ) = A > XF RAC ( ^ ) = 0 . 5475 

XFRACt 5 )=3 j XFRAC(6) = 2.0 
XFRAC ( 7 ) — 2 } XFRAC ( 8 ) =1 . A 



65 



GROUP 4. Y-DIRECTION GRIO SPECIFICATION 
NY = 12 



YFRACt l) = -3 


> 


YFRACt 2J=0.562 


YFRACt 3 ) = 2 


> 


YFRACt 4 )=0 . 314 


YFRAC ( 5 )=3 


) 


YFRAC t 6 ) =0 .562 


YFRAC( 7 )=4 


> 


YFRACt 8 )=0. 75 


GROUP 5. Z 
NZ=40 


-DIRECTION GRID SPECIFICATION 


ZFRACt 1)=-1 


> 


ZFRACt 2) = 1.0 


ZFRACt 3 ) = 1 


> 


ZFRACt 4 ) = 1 . 7 


ZFRACt 5 ) = 1 


> 


ZFRAC (6 ) = 2 . 0 


ZFRAC ( 7 )=1 


> 


ZFRAC t 8 )=1 . 7 


ZFRAC ( 9 ) = 2 


> 


ZFRACt 10 )=3.0 


ZFRAC( 11)=1 


\ 


ZFRACt 12 1 = 1.5 


ZFRACt 13 ) = 1 


> 


ZFRACt 141=1.0 


ZFRACt 15 ) = 1 


> 


ZFRACt 16 ) = 1 . 7 


ZFRACt 17 ) = 2 




ZFRACt 18 ) = 1 . 0 


ZFRACt 19 ) = 1 


> 


ZFRAC ( 20 ) = 1 . 7 


ZFRACt 21 ) = 2 


> 


ZFRACt 22 ) = 1 . 0 


ZFRACt 23 )=3 


> 


ZFRACt 24 )=0 .566666 


ZFRACt 25 )=2 


> 


ZFRACt 26 )=0 . 35 


ZFRACt 27 ) = 7 


> 


ZFRACt 28 )=3 . 0 


ZFRACt 29 )=1 


> 


ZFRACt 30 ) = 2 . 0 


ZFRACt 31 )=4 


) 


ZFRACt 32) =0.75 


ZFRACt 33 )=4 




ZFRACt 34 ) = 0. 575 


ZFRAC l 35 )=4 


> 


ZFRACt 36 )=0 .5 


ZFRACt 37 ) = 1 


> 


ZFRAC t 38 ) =3 . 0 



GROUP 6. BODY-FITTED COORDINATES OR GRID DISTORTION 

GROUP 7. VARIABLES STORED, SOLVED & NAMED 
SOLVE ( PI ,H1 >U1 >V1 >W1 ) 

SOLUTNI P1,Y,Y,Y,N,N,N) 

STORE ( TMP1,RH01 ) 

GROUP 8. TERMS (IN DIFFERENTIAL EQUATIONS) 2 DEVICES 

GROUP 9. PROPERTIES OF THE MEDIUM 
TAMB=295. 

TJET=1944. 

PJET=0. 

MJET=197.72 
PRESS0=101325 . 

RH01A=0. 

GASCON=286. 

RHOlBsl. /GASCON 

RHOJET=( ( PRESSO + PJET )*RH01B )/TJET 

WJET=MJET/RHOJET 

TMP1A=TINY 

*** TMP1B * 1. DIVIDED BY C-SUB-P) USED IN TMP1=GRND5 
CSU3P=1004. 

HAMB=TAMB*CSUBP 
TMP1B=1 ./CSUBP 
HJET=CSUBP*TJET 
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** NEXT 2 STATEMENTS INVOKE COMPRESSIBLE FLOW IN GREX1 
TMP1=GRND2 
RH01=GRN05 
DRH10P =GRN05 
ENUL-1 . E-05 
E NUT =200 . *ENUL 
TURMOO(KEMOOL) 

STORE (VIST) 

KEINJ=0.5*(0.005*HJET)**2 > EPINJ=0. 09*KEINJ**1. 5/1 0.05*0.5 ) 
KEINA=0.5*( .005*1.0 )**2 } EPINA=0. 09*KEINA**1 .5/( 0 . 1 ) 

GROUP 10. INTER-PHASE -TRANSFER PROCESSES AND PROPERTIES 

GROUP 11. INITIALIZATION OF VARIABLE OR POROSITY FIELDS 
FIINITI HI )=HAMB 

FIINIT ( TMP1 )=TAMB > FIINIT( DENI )=1 . 0 

FIINITI KE )=KEINA > FIINIT( EP ) = EPINA 

FIINIT( PI )=0 . 0 

FIINITI HI )=0.0 

FIINITI V1)=0.0 

FIINITI U1 ) = 0. 0 

PATCH! HINIT, INI VAL, 1,5, 1,8, 15, 30, 1,1) 

INI T( HINIT, HI ,0.0,0. 2*WJET ) 

I NIT I HINIT , HI ,0 . 0 ,CSUBP*0 . 2*TJET ) 

INITI HINIT ,KE , 0 . 0 ,KEINJ*0 . 2 ) 

INITI WINIT,EP,0.0 , EPINJ*0 . 2 ) 

*** SEE GROUP 15 FOR POROSITY FIELOS 

GROUP 12. CONVECTION ANO DIFFUSION ADJUSTMENTS 

GROUP 13. BOUNOARY CONDITIONS ANO SPECIAL SOURCES 

**OEFINE SOURCES AND SINKS FOR JET ENGINE** 

CONPORI 0.0, CELL, 1,1, 4,5, 12,14) 

PATCH I JE TIN , HIGH ,1,1,4,5,11,11,1,1 ) 

COVALI JETIN,P1,FIXFLU,-MJET) 

C OVA LI JETIN ,H1 ,ONLYM$ ,SAME ) 

PATCH! JETOUT, LOW, 1,1, 4, 5, 15,15,1,1) 

COVALI JETOUT, PI, FIXFLU,MJET) 

COVALI JETOUT, HI, ONLYMS,H JET) 

COVALI JETCUT , HI ,ONLYMS ,HJET ) 

COVALI JETOUT, KE,ONLYMS,KEINJ) 

COVALI JETOUT, EP,ONLYMS,EPINJ) 

**EAST BOUNOARY TO ATMOSPHERE** 

PATCH! OUTSIOE , EAST, 10, 10, 1,12,1, 40, 1, 1 ) 

COVALI OUTSIOE, PI, 0.1, 0.0) 

COVALI OUTSIOE ,H1 , FIXVAL ,HAMB ) 

**TOP OF HUSH HOUSE ANO TUBE INCLUDING TUBE LIP** 

CONPOR! 0.0, NORTH, 1,9, 12, 12, 1,19) 

CONPORI 0.0, NORTH, 1,5, -8, -8, 18, 2 7 ) 

CONPORI 0 . 0 , NORTH , 1 ,5 , -2 , -2 , 35,35 ) 

CCNPOR (0.0 , NORTH , 1 , 5 , -4 , -4 , 36 , 36 ) 

CONPORI 0.0, NORTH, 1,5, -5, -5, 37, 37 ) 

CONPORI 0.0, NORTH, 1,5, -6, -6, 38, 38) 
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CONPORl 0.0, NORTH, 1,5, -7, -7, 39, 39) 



**TOP OF MODEL TO ATMOSPHERE** 

PATCH ( TOPOLTT , NORTH , 1 , 10 , 12 , 12 , 20 ,40 , 1 , 1 ) 
COVALITOPOUT,P1,0.1,0.0 ) 

COVAL ( TOPOUT, Hi , ONLYMS ,HAMB ) 

PATCHlTOPP,NORTH,10,lO,12,12,l,19,l,l) 

COVAL(TOPP,P1,0.1,0.0) 

COVALITOPP, HI, ONLYMS, HAMB) 

** FLOOR OF HUSH HOUSE, TUBE AND GROUND** 

CONPORl 0.0 , SOUTH ,1,10 , 1 ,-1,1,40 ) 

**FORE AND AFT BULKHEADS IN HUSH HOUSE AND TUBE** 
CONPORl 0.0, LOW, 1,9, 1,12, -1,-1) 

CONPORl 0.0, LOW ,6 ,9,1,12,-20,-20) 

CONPORl 0 . 0 , LOW , 1 ,5 , 9 , 1 2 , -20 , -20 ) 

CONPORl 0 . 0 , LOW , 1 , 5 , 1 , 2 , -35 , -35 ) 

CONPORl 0.0, LOW ,1,5, 3,4, -3 6, -3 6) 

CONPOR( 0.0, LOW ,1,5, 5, 5, -3 7, -3 7) 

CONPORl 0.0, LOW , 1 ,5 ,6 ,6 , -38,-38 ) 

CONPOR( 0.0, LOW , 1 ,5 , 7 , 7 , -3 9 , -39 ) 

CONPORl 0.0, HIGH, 1,5, 8, 8, -3 9, -3 9) 

**VERY FRONT MOOEL STRIP TO ATMOSPHERE** 

PATCH( FRONT, LOW, 10, 10, 1,12, 1,1, 1,1) 

COVAL ( FRONT , PI, 0.1, 0.0 ) 

COVAL ( FRONT, HI, ONLYMS, HAMS ) 

**THE VERY ENO OF OUR MOOEL TO ATMOSPHERE** 

PATCH! ENOENO, HIGH, 1,10, 1,12, 40, 40,1,1) 

COVAL l ENDEND, PI, 0. 1, 0. 0 ) 

COVAL ( ENDEND , H 1 , ONLYMS , H AMB ) 

**0UTSI0E EAST WALL IN HUSH HOUSE** 

CONPORl 0 . 0 , E AST ,-9,-9,9,12,1,19) 

CONPORl 0.0, EAST, -9, -9,1, 8, 1, 1 ) 

CONPORl 0 . 0 , E AST ,-9,-9, 1,8, 3, 3) 

CONPORl 0.0, EAST, -9, -9, 1,8, 5, 8) 

CONPORl 0.0, EAST, -9, -9, 1,8, 10, 11) 

CONPORl 0.0, EAST ,-9,-9,1,8,13,14) 

CONPORl 0.0, EAST, -9,-9,1,8,18,19) 

**OUTSIDE OF TUBE PLUS TUBE LIP IN HUSH HOUSE** 
CONPORl 0.0, EAST, -5, -5,1, 8, 18, 39) 

GROUP 14. DOWNSTREAM PRESSURE FOR PARAB= . TRUE . 

GROUP 15. TERMINATION OF SWEEPS 
F$W£EP=1> LSWEEP=180 

GROUP 16. TERMINATION OF ITERATIONS 

GROUP 17. UNOER-RELAXATION DEVICES 
RE ALl OELT , FACT ) 
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DELT=50 . 0/1 N2*0 . 2*WJET ) 

FACT=10. 

FACT=2.5 
FACT -0 . 75 
FACT=3.5 
FACT=5. 

FACT=0 . 50 
FACT=0 . 25 

RELAX!P1,LINRLX,0.3) 

RE LAX l U1 , FALSDT , FACT*DE LT ) 

RE LAX! VI , FALSDT > FACT*DE LT ) 

RELAX! HI , FALSDT, FACT*DELT ) 

RE LAX! HI , FALSDT , 2 . 5*FACT*DE LT ) 

RELAX! KE, FALSDT, FACT*DELT) 

RELAX! EP, FALSDT, FACT*DELT) 

GROUP 18. LIMITS ON VARIABLES OR INCREMENTS TO THEM 
VARMIN! U1 )=-300 . 

VARMAX! U1 )=1000 . 

VARMIN! VI )=-300. 

VARMAX! VI ) = 1000. 

VARMIN! HI ) = -100 . 

VARMAX! HI) = 1500. 

VARMIN! TMP1)=250 
VARMAX! TMP1 )=TJET 
VARMIN! HI ) = 1000. 

VARMAX! HI )=CSUBP*TJET 
VARMIN! PI ) = -3 .5E+05 
VARMAX! PI )= 3.5E+05 
VARMIN! RH01)=. 05 
VARMAX! RHOl) =5. 

VARMIN! KE )=l.E-08 
VARMAX! KE )=1.3*KEJET 
VARMIN! EP ) S 1 . E -08 
VARMIN! VIST )=0.1 
VARMAX! EP )=1 . 3*! KEJET**1 .5 ) 

GROUP 19. DATA COMMUNICATED BY SATELLITE TO GROUND 

GROUP 20. PRELIMINARY PRINT -OUT 

GROUP 21. PRINT-OUT OF VARIABLES 

GROUP 22. SPOT-VALUE PRINT-OUT 

^MONITOR THE OUTPUT AT DESIGNATED POINTS** 
IXMON=l > IYM0N=4 i I2M0N=16 
TSTSWP=1 > LUPR3=6 
ABSI2=1 . 0 ) 0RSI2=1 . 0 

GROUP 23. FIELD PRINT-OUT AND PLOT CONTROL 
NUNCLS=10*NXPRIN 
ITABL=3 > NPLT =1 ) N2PRIN = 1 
IPLTL=LSWEEP 
NPRINT=3 
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GROUP 24. DUMPS FOR RESTARTS 
SAVE=T 
RESTRTt ALL) 

STOP 
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appendix b 

SELECTED PHOENICS OUTPUT 



W V W.M M.1L. 



»»»« » »«» » »« « 11 ******* ** 



CCCC HHH 
CCCCCCCC HHHHH 
CCCCCCC HHHHHHHHHH 
CCCCCCC HHHHHKHHHHHH 
CCCCCC HHHHHHHHHHHHH 
CCCCCCC HHHHHHHHHHHH 
CCCCCCC HHHHHHHHHH 
CCCCCCCC HHHHH 
CCCC HHH 



PHOENICS VERSION 1.3, 03 SEPT 1986 
1C) COPYRIGHT 1984 

CONCENTRATION HEAT AND MOMENTUM LTD 
ALL RIGHTS RESERVED. 

CHAM LTD, BAKERY HOUSE, 40 HIGH ST 
WIMBLEDON, LONDON, SW19 5AU 
TEL: 01-947-7651 > TELEX: 928517 
FACSIMILE: 01-879-3497 
THE OPTION LEVEL IS -18 



THIS CODE MAY ONLY BE USED UNDER THE TERMS AND CONDITIONS 
OF A LICENCE AGREEMENT WITH CHAM LTD. 

REPLICATION OF THIS CODE IS PROHIBITED UNLESS 



SPECIFICALLY AUTHORISED IN WRITING BY CHAM LTD. 



GREX1 OF 15/07/86 HAS BEEN CALLED 

FORMATTED SATLIT DATA READ FROM DF10 FOR IRUN= 1 

*** STORAGE INFORMATION «** 

F DIMNSN=200000 OCCUPIED=123707 ESTIMATED MINIMUM DIMNSN= 23267 
DEP VR8L= 67200 OLD VRSL = 0 3D COEFF= 23880 3D 0V0PS= 14400 

INITIAL FIELDS READ FROM DF09 

#** *** ********* * ** * ***** * * * * * ** *** ***** * * **** * *+H + ** *#-#***-** * 

GROUP 1. RUN TITLE & NUMBER 

3A V ^VV WW.YAW V W , -- 1Y ^VWWWWV VWW>/W I L/--<WWW WM^WWWW'UWVVA W W ^ VW w 



TEXT) MAS JACKSONVILLE HUSH CHAMBER ) 

M MJy. 1A V M x X W W V V W V 3A V %A W W W X! SA V 3A 3A W W W W W W 3A LY M M W W. W W U WWWWWWWWWW 

AA A AA A A A A A A A * A A A A A^T 

wvwwww ^.^v^v v wwv v wwwv www w w wwww ^ w w w w w w vr w ^ w w w w w w w w w w w WWW w w u w w w 
W AA^^AA A AAA A A^' AA AAA A A AA A^\ A A A A^3 A^> A AA A^TA A A A AA A^\A A AA AA A 

IRUNN = 1 

GROUP 2. TRANSIENCE 
STEADY = T 

GROUP 3. X-DIRECTION GRID SPACING 
CARTES = T 

NX = 10 

XULAST = 1 . OOOE+OO 

METHOD OF PAIRS USED FOR GRID SETTING 



XFRAC < 


1) 


= 


-1. OOOE+OO 


iXFRAC ( 


2 J 


= 


3 . 100E-01 


XFRAC ( 


3) 


= 


<* . OOOE+OO 


)XFRAC ( 


4 ) 


= 


5.475E-01 


XFRAC i 


5) 


s 


3. 000E+00 


jXFRAC ( 


6 ) 


= 


2. 000E+00 


XFRAC ( 


7) 


= 


2. 000E+00 


>XFRAC ( 


8 ) 


= 


1.400E+00 



*»*»*»»»*»«* * *»*»**»tt^»»*****»»*^** *»**»»* 

GROUP <♦. Y-DIRECTION GRID SPACING 
NY =12 

YVLAST = 1.000E+00 

METHOD OF PAIRS USED FOR GRID SETTING 
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YFRAC 


( 


1) = 


-3 . 000E *00 


jYFRAC ( 


2) 


s 


5. 620E-01 


YFRAC 


( 


3) = 


2.000E+00 


jYFRAC ( 


4) 


= 


3 . 140E-01 


YFRAC 


( 


5) = 


3.000E+00 


jYFRAC ( 


6) 


= 


5. 620E-01 


YFRAC 


( 


7) = 


4.000E+00 


jYFRAC ( 


8) 


= 


7.500E-01 



2 - 



GROUP 5. 
PARAB = 

N2 = 

ZWLAST = 1 



DIRECTION GRID SPACING 
F 

40 

000E+00 



METHOD OF 


PAIRS USED FOR GRID 


SETTING 




2FRAC 


( 


1) 


= 


-1 . 000E ♦00 


J2FRAC 


( 


2) 


s 


1 . 000E ♦OO 


2FRAC 


( 


3) 


= 


1.000E+00 


J2FRAC 


( 


4) 


s 


1 . 700E 4-00 


2FRAC 


( 


5) 


= 


1.000E+00 


J2FRAC 


( 


6) 


= 


2.000E+00 


2FRAC 


( 


7) 


s 


1.000E+00 


I2FRAC 


( 


8) 


= 


1.700E+00 


2FRAC 


( 


9) 


= 


2.000E+00 


&2FRAC 


( 


10) 




3.000E+00 


2FRAC 


( 


11) 


= 


1 . 000E *00 


J2FRAC 


( 


12) 


= 


1.500E+00 


2FRAC 


( 


13) 


= 


1.000E+00 


>2FRAC 


( 


14) 




1.000E+00 


2FRAC 


( 


15) 


= 


1.000E+00 


J2FRAC 


( 


16) 


= 


1. 700E+00 


2FRAC 


( 


17) 


= 


2 . 000E *00 


J2FRAC 


( 


18) 


= 


1 . 000E ♦CO 


2FRAC 


( 


19) 


= 


1.000E+00 


>2FRAC 


( 


20) 


= 


1 . 700E+00 


2FRAC 


( 


21) 


= 


2.000E+00 


J2FRAC 


( 


22) 


= 


l.OOOEfOO 


2FRAC 


( 


23) 


= 


3.000E+00 


)2FRAC 


( 


24) 


= 


5.667E-01 


2FRAC 


( 


25) 


= 


2 . OCOE + 00 


*2FRAC 


( 


26) 


= 


3.500E-01 


2FRAC 


( 


27) 


= 


7.000E+00 


J2FRAC 


( 


28) 


= 


3.000E+00 


2FRAC 


( 


29) 


= 


1 . 000E +00 


J2FRAC 


( 


30) 


= 


2 . 000E *00 


2FRAC 


( 


31) 


= 


4 . OOOE 4-00 


>2FRAC 


( 


32) 


= 


7.500E-01 


2FRAC 


( 


33) 


= 


4.000E+00 


>2FRAC 


( 


34) 


= 


5. 750E-01 


2FRAC 


( 


35) 


= 


4.000E+00 


J2FRAC 


( 


36) 




5.000E-01 


2FRAC 


( 


37) 




l.OOOEtOO 


J2FRAC 


( 


38) 


= 


3 . 000E 4-00 



GROUP 6. BFC AND GEOMETRY CHANGES 
«****-***#*-**#*« 

GROUP 7. VARIABLES: STORE 0 , SOLVE 0 > NAMED 
ONEPHS = T 

NAME ( 1) =P1 jNAME ( 3) =U1 
NAME ( 5) -VI )NAME( 7) =H1 

NAME (12) -KE jNAME(13> =EP 

NAME ( 14 ) sHl jNAME ( 19 ) -HPOR 

NAME ( 20 ) =NPCR ;NANE(21) =EPOR 
NAME (22) =VPOR »f4AME (23) =VIST 
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NAME (24) = RH01 >NAME<25) =TMP1 
SOLlTTNfPl ,Y,Y,Y,N,N,N) 

SOLUTNtUl ,Y,Y,N,N,N,N) 

SOLUTN( VI ,Y,Y,N,N,N,N) 

SOLUTNt HI ,Y,Y,N,N,N,N) 

SOLUTN(KE ,Y,Y,N,N,N,N) 

SOLUTN(EP ,Y,Y,N,N,N,N) 

SOLUTN( HI ,Y,Y,N,N,N,N) 
SOLUTN(HPOR,Y,N,N,N,N,N) 
SOLUTN(NPOR,Y,N,N,N,N,N) 
SOLUTN(EPOR,Y,N,N,N,N,N) 

SOLUTN( VPOR,Y,N,N,N,N,N) 

SOLUTW VIST,Y,N,N,N,N,N) 

SOLUTNl RH01,Y>N,N,N,N,N) 
SOLUTN(TMPl,Y,N,N,N,N,N) 

DENI = 24 

VIST = 23 

EPOR s 21 jHPOR = 19 >NPOR 

TEMPI = 25 



20 >VPOR 



22 



GROUP 8. TERMS * DEVICES 
TERMS (PI ,Y,Y,Y,N,Y,N) 
TERMS (U1 ,Y,Y,Y,N,Y,N) 

TERMS (VI ,Y,Y,Y,N,Y,N) 
TERMS (HI ,Y,Y,Y,N,Y,N) 
TERMS (KE ,N,Y,Y,N,Y,N) 
TERMS (EP ,N,Y,Y,N,Y,N) 
TERMS (HI ,Y,Y,Y,N,Y,N) 
DIFCUT = 5 . Q00E-01 >ZDIFAC 
GALA * F JAODOIF = 

NEWRH1 = T 

NEHENT = T 



= 1.000E+00 
F 



GROUP 9. PROPERTIES 



RHOl 


S* 


-1 . 016E + 04 


jTMPl 


=-1.013E+04 


jELl 


s-1.01SE + 0<+ 


RH01A 


= 


O.OOOE+OO 


>RH01B 


* 3.496E-03 


JRH01C 


= 7 . 143E-01 


PRESSO 

TMP1A 


= 


1.013E+05 

1.000E-20 


>TMP1B 


= 9. 960E-Q4 


jTMPIC 


= O.OOOE+OO 


TEMPO 

ELIA 


= 


0.000E+00 

O.OOOE+OO 


,EL1B 


= O.OOOE+OO 


> EL1C 


= O.OOOE+OO 
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ENUL = 1.000E 


-05 jENUT 


=-1.014E*04 






ENUTA = O.OOOE+OO >ENUTB 


= 0 . 000E+00 


>ENUTC * O.i 


0RH1DP =-1.016E*04 










PRNDTLtUl ) = 


1.000E+00 


jPRNDTLt VI 


) 


= 


1 . 000E *00 


PRNDTLl HI ) = 


1 . 000E+00 


>PRNDTL( KE 


) 


= 


1 . 000E *00 


PRN0TKEP ) = 


1 . 000E +00 


jPRNDTLt HI 


) 


= 


1 . 000E + 00 


PRT (U1 J = 


1.000E+00 


>PRT tVl 


) 


= 


1 . 000E *00 


PRT (HI ) = 


1 . 00CE+00 


)PRT (KE 


) 


s 


1 . OOOE +00 


PRT tEP ) = 


1 . 314E *00 


>PRT (HI 


) 


= 


1 . 000E ♦OO 



OOOE+OO 



htWWUW V WV V W W W W W W W ^ W V/ V V W WV WW W WV^V W WWwV^f X V U U U WWVhfVVt. ~W • I WWW^MWW U 

n n inc “ w w jrinrww k x ttit k f w kw k wtc h w m irwinc w wn ^ a n ^ hr !T wk k k x h k w 

GROUP 10. INTER-PHASE TRANSFER PROCESSES 



GROUP 11. INITIALIZE VAR/POROSITY FIELOS 



FI INI T( PI ) = 
FIINIT( VI ) = 
FIINIT ( KE ) = 
FIINITC HI ) = 
FIINITt NPOR J = 
FIINITt VPOR ) = 
FIINITt RH01) s 



■1 . 022E *04 
■1 . 022E *04 
•1.022E+04 
■1. 022E *04 
*1 . 022E *04 
•1 . 022E +04 
*1. 022E+04 



> FIINITt U1 ) 
>FIINITtWl J 
>FIINIT(EP ) 
jFIINIT(HPOR) 
>FIINIT( EP0R ) 
)FIINITC VIST J 

> FIINITt TMP1 ) 



-1 . 022E ^04 
-1 . 022E *04 
-1 . 022E *04 
-1.022E+04 
-1 . 022E ^04 
-1.022E+04 
-1.022E+04 



PATCHt S0UTH1 ,INIVAL, 1, 1, 3, 



INlTt SOUTHl 


,V1 , 


0. 000E+00 , 


0.000E + 00 ) 


INIT( SOUTHl 


,NPOR , 


0 . 000E + 00 , 


0 . 000E + 00 > 


PATCHt L0H1 


,INI VAL , 1, 


1, 4, 5 


INITtLOWl 


,H1 , 


0 . OOOE *00 , 


0 . OOOE ♦OO J 


INITt L0KI 


, HPOR , 


0 .000E+00 , 


0 . OOOE tOO ) 


PATCHt CELL1 


,INIVAL, 1, 


1, 4, 5 


INITt CELL1 


,P1 , 


0 . OOOE ♦OO , 


0 . 000E + 00 J 


INlTt CELL1 


,U1 , 


0 . OOOE ♦OO , 


0.000E+00 ) 


INITt CELL1 


>V1 , 


0 . OOOE *00 , 


0 . OOOE *00 ) 


INITtCELLl 


>W1 , 


0 . OOOE +00 , 


0 . 000E + 00 ) 


INITt CELL1 


,KE , 


0 . OOOE *00 , 


0 . OOOE +00 ) 


INITtCELLl 


,EP , 


O.OOOE+OO, 


0 . OOOE ♦OO ) 


INITtCELLl 


,H1 , 


0 .OOOE '►OO , 


0 . OOOE ♦OO ) 


INITtCELLl 


,HPOR, 


0 . OOOE ♦OO , 


0 . OOOE + 00 ) 


INITtCELLl 


,NP0R, 


0 . OOOE +00 , 


0 . OOOE +00 J 


INITtCELLl 


,EP0R , 


0 . OOOE+OO , 


0 . OOOE +C0 ) 



12 , 14 , 



11 , 11 , 



12, 14, 



1 ) 



1 ) 



1 ) 
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PATCH) EAST21 ,INIVAL, 9, 9, 1, 8, 

INITI EAST21 ,U1 , O.OOOE+OO, O.OOOE+OO) 

INIT! EAST21 ,EPOR, O.OOOE+OO, O.OOOE+OO) 


3, 


3, 


1, 


1) 


PATCH! EAST22 ,INIVAL, 9, 9, 1, 8, 

INITI EAST22 ,U1 , O.OOOE+OO, O.OOOE+OO) 

INITI EAST22 ,EPOR, O.OOOE+OO, O.OOOE+OO) 


s. 


8, 


1, 


1) 


PATCH! EAST23 ,INIVAL, 9, 9, 1, 8, 

INITI EAST23 ,U1 , O.OOOE+OO, O.OOOE+OO) 

INITI EAST23 ,EPOR, O.OOOE+OO, O.OOOE+OO) 


10, 


11, 


1, 


1) 


PATCH! EAST24 ,INIVAL, 9, 9, 1, 8, 

INIT(EAST2<+ ,U1 , O.OOOE+OO, O.OOOE+OO) 

INITI EAST2<+ ,EPOR, O.OOOE+OO, O.OOOE+OO) 


13, 


14, 


1, 


1) 


PATCH! EAST25 ,INIVAL, 9, 9, 1, 8, 

INITI EAST25 ,U1 , O.OOOE+OO, O.OOOE+OO) 

INITI EAST25 ,EPOR, O.OOOE+OO, O.OOOE+OO) 


18, 


19, 


1, 


1) 


PATCH! EAST26 ,INIVAL, 5, 5, 1, 8, 

INITI EAST26 ,U1 , O.OOOE+OO, O.OOOE+OO) 

INIT! EAST26 ,EPOR, O.OOOE+OO, O.OOOE+OO) 


18, 


39, 


1, 


1) 


INIAOD * F 

LUFI = 9 

NAMFI *DF09 

W yg .w w W V W ^ V/ W W \/ w w W W W W X ^ M w W W W W UWWWWWWVW w ^ X w w W 

WKW F W W 1r1rlr^KK7^ A R n R K * K r K Kfl R R R R R R R TTR R R R WWW K R R IT 








• 


GROUP 13. BOUNDARY 2 SPECIAL SOURCES 










PATCH ( KESOURCE ,PHASEM, 1, 10, 1, 12, 

COVAU KESOURCE ,KE , -1 . 015E *04, -1 . OISE *04 ) 
COVAU KESOURCE ,EP ,-l .015E+04,-l .015E+04 ) 


1, 


40, 


1, 


1) 


PATCH! JETIN ,HIGH , 1, 1, 4, 5, 

C0VAL( JETIN ,P1 , 1.000E-10,-1.977E*02) 

COVAUJETIN ,H1 , 0 . 000E *00 , -1.026E+04 ) 


11, 


11, 


1, 


1) 


PATCH( JETOUT ,LOW , 1, 1, 4, 5, 

COVAUJETOUT ,P1 , 1.000E-1G, 1.977E+02) 


IS, 


IS, 


1, 


1) 
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COVAL! HWALL12 


,U1 


> 


COVAL! HWALL12 


,V1 


9 


COVAL! HWALL12 


#KE 


9 


COVALt HHALL12 


»EP 


9 


PATCH! LKALL12 


,LHALL 


COVAH LKALL12 


,U1 


9 


COVALl LHALL12 


,V1 


9 


COVAL! LWALL12 


,KE 


9 


COVAH LHALL12 


,EP 


9 


PATCH! HHALL13 


, HWALL 


COVAL! HWALL13 


,U1 


9 


COVAL! HWALL13 


,vi 


9 


COVAL! HHALL13 


>KE 


9 


COVAL! HHALL13 


»EP 


9 


PATCH! LHALL13 


,LWALL 


COVAH LWALL13 


»U1 


9 


COVAH LWALL13 


,V1 


9 


COVAH LHALL13 


,KE 


9 


COVAL! LWALL13 


,EP 


9 


PATCH! HWALL14 


, HWALL 


COVAH HWALL14 


,U1 


9 


COVAL! HWALL 14 


,V1 


9 


COVAH HUALL14 


,KE 


9 


COVAL! HWALL14 


,EP 


9 


PATCH! LWALL14 


,LWALL 


COVAL! LWALL14 


»U1 


9 


COVAH LWALL14 


,V1 


9 


COVAH LWALL14 


,KE 


9 


COVAL! LKALL14 


,EP 


9 


PATCH! HWALL15 


,HWALL 


COVAL! HWALL 15 


fUl 


9 


COVAH HWALL15 


,V1 


9 


COVAL! HWALL1S 


iKE 


9 


COVAL! HHALL 15 


,EP 


9 



-1.013Et04, O.OOOEtOO) 
-1 . 013E *04 , 0 . OOOE +00 ) 
-1.013Et04,-1.0l3Et04) 
-1.013Et04,-1.013Et04) 

, 1, 5, 9, 12, 
-1.013E+04, 0 . OOOE *00 ) 
-1 . 013E +04 , 0. OOOE tOO ) 
-1.013Et04,-1.013Et04 J 
-1.013Et04,-1.013Et04) 

, 1, 5, 1, 2, 
-1 . 013E t04 , 0 . OOOE tOO ) 
-1.013Et04, 0 .OOOE +00 ) 
-1.013Et04,-1.013Et04 ) 
-1.013Et04,-1.013Et04) 

, 1, 5, 1, 2, 
-1 . 013E +04 , 0 . OOOE + 00 ) 
-1 . 013E t04 , 0 . OOOE tOO ) 
-1.013Et04,-1.013Et04 ) 
-1.0l3Et04,-1.013Et04) 

» 1, 5 » 3, 4 , 
-1 . 013E t04 , 0 . OOOE + 00 ) 
-1.013Et04, O.OOOEtOO) 
-1.013Et04,-1.013Et04 ) 
-1.013Et04,-1.013Et04) 

, 1, 5, 3, 

-1 .013E+04 , O.OOOEtOO ) 
-1 . 013Et04, O.OOOEtOO) 
-1.013Et04,-1.013Et04 ) 
-1.013Et04,-1.013Et04) 

> 1> 5, 5, 5, 
-1 . 013E t04 > O.OOOEtOO) 
-1.013EtO<t, O.OOOEtOO) 
-1.013Et04,-1.013E+04) 
-1.0l3Et04,-l.Cl5Et04 ) 



20 , 20 , 1 , 1 ) 



34, 34, 1, 1) 



35, 35, 1, 1) 



35, 35, 1, 1) 



36, 36, 1, 1 ) 



36, 36, 1, 1) 
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PATCH! LWALL15 


, LHALL , 1, 5, 5, 5, 


37, 


37, 


If 


1) 


COVAL! LHALL15 


,U1 


,-1.013E+04, 0 . 000E *00 ) 










COVAL! LKALL15 


,V1 


»-1.013E+0 < *» O.OOOE + OO) 










COVAL! LHALL15 


,KE 


,-1.013E+04,-1.013E*04) 










COVAL! LHALL15 


>EP 


,-1.013E+04,-1.013E+04) 










PATCH! HWALL16 


,HWALL , 1, 5, 6, 6, 


37, 


37, 


l, 


1) 


COVAL! HHALL16 


,U1 


,-1.013E404, o.oooe*ooj 










COVAL! HHALL16 


,V1 


,-1.013Et04, 0 • 000E *00 ) 










COVAL! HWALL16 


,KE 


,-1.013E*04,-1.013E+04) 










COVAL! HKALL16 


,EP 


,-1.013E*04,-1.013E*04) 










PATCH! LWALL16 


> LHALL , 1, 5, 6, 6, 


38, 


38, 


I, 


1 ) 


COVAL! LWALL16 




,-1.013E*04, 0.000E+00) 










COVAL! LWALL16 


,vi 


,-1.013E*04, 0.000E+0Q) 










COVAL! LHALL16 


>KE 


,-1.013E*04,-1.013E*04) 










COVAL! LKALL16 


,EP 


,-1.013E*04,-1.013E*04) 










PATCH! HWALL17 


>HWALL , 1, 5, 7, 7, 


38, 


38, 


If 


1) 


COVAL! HHALL17 


,U1 


,-1.013E*04, O.OOOE+OO) 










COVAL! HWALL17 


,V1 


, -1 . 013E *04 , O.OOOEtOOJ 










COVAL! HHALL17 


>KE 


#•1. 013E*04,-1.013E*04 ) 










COVAL! HHALL17 


>EP 


,-1.013E*04,-1.013E*04) 










PATCH! LKALL17 


, LHALL , 1, 5, 7, 7, 


39, 


39, 


If 


1) 


COVAL! LHALL17 


>U1 


9 -1 . 013E ^04 9 O.OOOE+OO) 










COVAL! LWALL17 


,V1 


,-1.013E*04, 0 . OOOE + OO ) 










COVAL! LWALL17 


>KE 


,-1.013E*04,-1.013E*04) 










COVAL! LKALL17 


>EP 


,-1.013E*04,-1.013E*04) 










PATCH! HWALL18 


,HHALL , 1, 5, 8, 8, 


39, 


39, 


If 


1) 


COVAL! HWALL13 


,U1 


,-1.013E*04, O.OOOEfOO) 










COVAL! HWALL18 


,V1 


,-1.013E*04, O.OOOE+OO) 










COVAL! HWALL18 


,KE 


,-1.013E*04,-1.013E*04) 










COVAL! HWALL18 


>EP 


,-1.013E*04,-1.013E*04) 










PATCH! LWALL18 


> LHALL , 1, S, 8, S, 


40, 


40, 


If 


1) 


COVAL! LHALL18 


,U1 


>-l . 013E ♦0<+ 9 0 . OOOE *00 ) 










COVAL! LHALL18 


,V1 


>-l . 013E+04 9 O.OOOEtOOJ 
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COVAL( LWAL118 


,KE 


> • 


COVAU LHALL18 


>EP 


i “ 


PATCH! FRONT 


,LOW 




COVAU FRONT 


,P1 


9 


COVAU FRONT 


> HI 


9 


PATCH! ENDEND 


,HIGH 


COVAU ENDEND 


,P1 


9 


COVAU ENDEND 


>H1 


9 


PATCH! EWALL19 


, EWALL 


COVAU EWALL19 


,V1 


9 ~ 


COVAU EHALL19 


,W1 


9 ” 


COVAU EWALL19 


,KE 


9 “ 


COVAU EWALL19 


,EP 


9 ~ 


PATCH! KWALL19 


,NWALL 


COVAL! KWALL19 


,V1 


9 ~ 


COVAU HHALL19 


,W1 


9 “ 


COVAU WWALL19 


,KE 


9 ~ 


COVAU WWALL19 


,EP 


9 " 


PATCH! EWALL20 


, EWALL 


COVAU EKALL20 


,V1 


9 * 


COVAL! EWALL20 


,H1 


9 ~ 


COVAL! EHALL20 


,KE 


9 * 


COVAL! EWALL20 


,EP 


9 * 


PATCH! WWALL20 


,WWALL 


COVAL! WWALL20 


,V1 


9 ~ 


COVAL! HHALL20 


,W1 


9 " 


COVAU WWALL20 


,KE 


9 ~ 


COVAU WHALL20 


» EP 


9 * 


PATCH! EHALL21 


, EWALL 


COVAL! EWALL21 


,V1 


, * 


COVAU EWALL21 


,W1 


9 ” 


COVAL! EHALL21 


,KE 


9 * 


COVAL! ENALL21 


> EP 


9 " 



1.013E+04,-1.013E*04) 

1.013E*04,-1.013E*04) 

, 10, 10, 1, 12, 

1.000E-01, 0 . 000E +00 ) 
0 . OOOE+OO , 2 . 962E+05 ) 

, 1, 10, 1, 12, 

1.000E-01, 0 . OOOE ♦OO ) 

0 . OOOE+QO , 2 . 962E+05 ) 

, 9, 9, 9, 12, 

1.013E+04, 0 . OOOE+OO ) 

1 . 013E+04 , O.OOOE+OO) 
1.013E*04,-1.013E*04) 
1.013E*04,-1.013E+04) 

, 10, 10, 9, 12, 

1 . 013E +04 , 0 . OOOE+OO ) 
1 . 0 13E ♦04 , 0 . OOOE ♦OC ) 
1.013E*04,-1.013E*04) 
1.013E*04,-1.013E*04) 

, 9, 9, 1, 8, 

1 . 013E+04 , O.OOOE+OO) 
1 . 013E *04 , O.OOOE+OO) 
1.013E*04,-1.013E*04 ) 
1.013E*04,-1.013E*04) 

, 10, 10, 1, 8, 

1 . 013E +04 , 0 . OOOE 4-00 ) 
1.015E+04, O.OOOE+OO ) 
1.013E*04,-1.013E*04) 
1.013E*04,-1.013E*04) 

, 9, 9, 1, 8, 

1 . 013E *04 , 0 . OOOE *00 ) 
1 . 013E +04 , 0 . OOOE *00 ) 
1 . 013E ♦04 ,-1.013E*04 I 
1.013E*04,-1.013E*0 4) 



1 , 1 , 1 , 1 ) 

40 , 40 , 1 , 1 ) 

1 , 19 , 1 , 1 ) 

1 , 19 , 1 , 1 ) 

1 , 1 , 1 , 1 ) 

1 , 1 , 1 , 1 ) 

3 , 3 , 1 , 1 ) 
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PATCH! WWALL21 


,HWALL 


COVAK KHALL21 


>V1 


> 


COVAU WWALL21 


>W1 


> 


COVAL! WWALL21 


>KE 


y 


COVAU WWALL21 


,EP 


9 


PATCH! EHALL22 


,EWALL 


COVAU EWALL22 


,V1 


9 


COVAU EWALL22 


>W1 


9 


COVAL! EWALL22 


>KE 


9 


COVAU EWAL122 


,EP 


9 


PATCH! HHALL22 


,HWALL 


COVAU WWALL22 


,V1 


9 


COVAU WWALL22 


,W1 


9 


COVAL! WWALL22 


,KE 


9 


COVAU KHALL22 


>EP 


9 


PATCH! EWALL23 


,EWALL 


COVAL! EHALL23 


,V1 


9 


COVAU EHALL23 


>W1 


9 


COVAU EHALL23 


>KE 


9 


COVAL! EWALL23 


,EP 


9 


PATCH! WWALL23 


,WWALL 


COVAU WHALL23 


>V1 


9 


COVAU WWALL23 


,H1 


9 


COVAU HHALL23 


,KE 


9 


COVAU WWALL23 


,EP 


9 


PATCH! EWALL24 


,EHALL 


COVAU EWALL24 


>V1 


9 


COVAL! EHALL24 


>W1 


9 


COVAU EWALL24 


,KE 


9 


COVAL! EKALL24 


>EP 


9 


PATCH! KWALL24 


,WWALL 


COVAL! WWALL24 


,V1 


9 


COVAL! !*'WALL24 


,W1 


9 



> 10 , 10 , 1 y 8 , 

-1.013E+04, 0 . 000E+00 ) 
-1.013E+04, 0 ♦000E+00 ) 
-1.013E*04,-1.013E*04) 
-1.013E*04,-1.013E*04) 

> 9, 9, 1, 8, 

-1.013E+04, 0.000E+00) 
-1 . 013E *04, 0 . 000E+00 ) 
-1.013E*04,-1.013E*04) 
-1 . 013E+04, -1 . 013E +04 ) 

, 10 , 10 , 1 , 8 , 

-1.013E+04, 0 . 000E *00 ) 
-1.013E+04, O.OOOE+OO) 
-1.013E*04,-1.013E*04) 
-1.013E*04,-1.013E*04) 

, 9, 9, 1, 8, 

-1.013E+04, O.OOOE+OO) 
-1.013E+04, 0 .OOOE+OO ) 
-1 . 013E *04,-1 . 013E +04 ) 
-1.013E*04,-1.013E*04) 

, 10 , 10 , 1 , 8 , 

-1 . 013E *04 , O.OOOEtOO) 
-1.0132*04, O.OOOE*OOJ 
-1.013E*04,-1.013E*04) 
-1.013E*04,-1.013E*G4) 

, 9, 9, 1, 8, 

-1.013E+04, 0.000E*00) 
-1.013E*04, 0 . OOOE tOO ) 
-1.013E*04,-1.013E*04) 
-1. 0 13 E* 04,-1 . 013E*04 ) 

, 10, 10, 1, 8, 

-1.013E+04, 0 . OOOE *00 ) 
-1.013E+04, 0 .000E*00 ) 



3, 3, 1, 1) 



5, 8, 1, 1) 



5, 8, 1, 1J 



10 , 11 , 1 , 1 ) 



10 , 11 , 1 , 1 ) 



13, 14, 1, 1) 



13, 14, 1, 1) 
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COVAU WWALL24 ,KE 
C0VAUWWALL24 ,EP 



, -1 . 013E*04,-1.013E*04 ) 
,-1.013E*04,-1.013E*04) 



PATCH( EWALL25 
COVAU EWALL25 
COVAU EHALL25 
COVAU EWALL25 
COVAU EWALL25 



,EHALL , 9, 9, 1, 8, 
,V1 ,-1.013E*04, 0 . 000E+00 ) 
,W1 ,-1.013E*04, 0.000E+00) 
,KE ,-1.013E*04,-1.013E*04) 
,EP ,-1.013E*04,-1.013E*04) 



18, 



19, 



1 , 



1) 



PATCH(WWALL25 
COVAU HWALL25 
COVAU WWALL25 
COVAU WWALL25 
COVALl WHALL25 



,KWALL , 10, 10, 1, 8, 

,V1 ,-1.013E*04, 0 .OOOE >00 ) 

,W1 ,-1.013E*04, 0 . 000E +00 ) 

,KE , -1 . 013E *04 , -1 . 013E *04 ) 
,EP ,-1.013E*04,-1.013E*04) 



18, 



19, 



1 , 1 ) 



PATCH! EWALL26 
COVAU EWALL26 
COVAU EWALL26 
COVALl EWALL26 
COVAU EWALL26 



,EWALL , 5, 5, 1, 8, 
,V1 ,-1.013E*04, 0.000E+00) 
,W1 ,-l .013E+04 , 0 .000E *00 ) 
,KE ,-1.013E*04,-1.013E*04) 
,EP , -1. 01 3E* 04,-1. 013E *04 ) 



18, 59, 



1, 



1) 



PATCHl WWALL26 
COVAU WWALL26 
COVAU WWALL26 
COVAU HWALL26 
COVAU HHALL26 
XCYCLE = 
KELIN = 



,KWALl , 6, 6, 1, 8, 

,V1 , -1 . 013E *04 , 0 . 000E *00 ) 

,H1 ,-1.013E*04, 0 . OOOE +00 ) 

,KE , -1.0 13 E* 04 , -1.013E + Q4 ) 
,EP ,-1.013E*04, -1. 013E +04 ) 
F 

0 



18, 39, 



1 , 



1) 



GROUP 14. DOWNSTREAM PRESSURE FOR PARAB 






GROUP IS. 
LSWEEP = 
LITHYD = 
RESREFl PI 
RESREFl VI 
RESREFl KE 
RESREFl HI 



TERMINATE SWEEPS 
180 >ISWC1 = 1 

1 >LITFLX = 1 > LITC 

) = 1 . OOOE -06 jRESREFl U1 ) 

) = 1 . 000E-06 >RESREFIW1 J 

) = 1 . OO0E-O6 jRESREFl EP J 

) = 1 . 000E-06 



1 >ITHC1 = 
1 . OOOE -06 
1. OOOE -06 
1 . OOOE -06 



GROUP 16. TERMINATE ITERATIONS 



1 
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LITER 


1 PI 


) 


= 


20 jLITER (U1 


) * 


20 




LITER 


(VI 


) 


= 


20 jLITER (HI 


) = 


20 




LITER 


(KE 


) 


S 


20 jLITER (EP 


) * 


20 




LITER 


! HI 


) 


= 


20 








ENDIT 


(PI 


) 


s 


1 . OOOE -06 JENDIT 


(U1 


) = 


1 . 000E-06 


ENDIT 


(VI 


) 


= 


1.000E-06 JENDIT 


(W1 


) = 


1 . 000E-06 


ENDIT 


(KE 


) 


= 


1.000E-06 JENDIT 


(EP 


) = 


1 . OOOE -06 


ENDIT 


(HI 


) 


= 


1 . 000E-06 









GROUP 17. 
RELAX! PI 
RELAX! U1 
RELAX! VI 
RELAX! HI 
RELAXtKE 
RELAX! EP 
RELAX! HI 
RELAX! VIST 
RELAX! RHOl 
RELAX! TNP1 



RELAXATION 
,LINRLX, 3.000E-01) 
, FALSDT , <► . 321E-03 ) 
, FALSDT , 4.321E-03) 
, FALSDT, 4 . 321E-03 ) 
, FALSDT, 4.321E-03) 
, FALSDT, <► . 321E-03 ) 
,FALSDT , 1.080E-02) 
,FALSDT , 1.000E+00) 
,FALSDT , 1 . OOOE+OO ) 
,FALSDT , 1 . OOOE+OO ) 



GROUP 18. 


LINITS 


VARNAX! PI 


) 


s 


3.500E+05 


VARNAX! U1 


) 


= 


1 . 000E+03 


VARNAX! VI 


) 


= 


1.000E+03 


VARNAX! HI 


) 


= 


1.500E+03 


VARNAX! KE 


) 


= 


1.000E+10 


VARMAXIEP 


) 


= 


1.000E+10 


VARNAX! HI 


) 


= 


1 . 952E *06 


VARNAX! HPOR) 


= 


1 . 000E+10 


VARNAX! NPOR) 




1 . 000E+10 


VARNAX! EP0R) 


= 


1 . 000E+10 


VARMAX! VPOR) 


= 


1.000EU0 


VARNAX! VIST) 




1.000E+10 


VARNAX! RHOl) 


= 


5 .000E + 00 


VARNAX! TNP1) 


= 


1 . 9^E*03 



jVARNIN! PI ) =-3.500E*05 
JVARMIN! U1 ) =-3.000E*02 
jVARNIN! VI ) --3 . OOOE *02 
JVARNIN! HI ) =-1.000E*02 
JVARNIN! KJE ) - 1.000E-10 
jVARNIN! EP ) = 1.000E-10 
jVARNIN! HI ) = 1.000E+03 
JVARMIN! HPOR ) =-1.000E*10 
jVARNIN! NPOR ) =-1.000E*10 
JVARNIN! EPORJ =-1.000E*10 
JVARNIN! VPOR) s-l.OOOEUO 
jVARNIN! VIST ) • 1.000E-20 
jVARNIN! RHOl) = 5.000E-02 
jVARNIN! TMP1 ) = 2.500E+02 






w v/ v w WW I/JW AZ I/WVW VU W WW W WwWWWWWWV^ w wwwwwwwwi_ w-_ W > r W W ^ W W W w • W W w W W V u M w 

* ^ ™ w w * w w w * w * 



GROUP 19. EARTH CALLS TO GROUND STATION 
USEGRD = T jUSEGRX - T 
NAMGRO =N0NE 
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GENK = T 



GROUP 20. PRELIMINARY PRINTOUT 
ECHO = T 



GROUP 21. 
INIFLD = 
OUTPUT ( PI 
OUTPUT! U1 
OUTPUT! VI 
OUTPUT! HI 
OUTPUT! KE 
OUTPUT! EP 
OUTPUT! HI 



PRINT-OUT OF VARIABLES 



F jSUBKGR = 
,Y,N,N,Y,Y,Y) 
,Y,N,N,Y,Y,Y) 
,Y,N,N,Y,Y,Y) 
,Y,N,N,Y,Y,Y) 
,Y,N,N,Y,Y,YJ 
,Y,N,N,Y,Y,Y) 
,Y,N,N,Y,Y,Y) 
OUTPUT! HPOR,N,N,N,N,N,N) 
OUTPUT! NPOR,N,N»N,N,N,N) 
OUTPUT! EPOR,N,N,N>N,N,N) 
OUTPUT! VPOR,N,N,N,N,N,N) 
OUTPUT! VIST, Y,N,N,N,N,N) 
OUTPUT! RH01,Y,N,N,N,N,N) 
OUTPUT! TMP1,Y>N,N,N,N,N) 



GROUP 22. MONITOR PRINT -OUT 
IXMON = 1 jIYKON = A >IZMON = 16 

NPRMON = 10000 jNPRMNT = 10000 >TSTSWP = 1 

HIGHLO = F 

GROUP 25 . FIELD PRINT-OUT i PLOT CONTROL 
NPRINT = 10000 iNUMCLS = 10 

NXPRIN = 1 > IXPRF = 1 >IXPRL = 10 

NYPRIN = 1 jIYPRF = 1 jIYPRL = 12 

NZPRIN = 1 jIZPRF = 1 jIZPRL = AO 

X2PR = F iYZPR = F 

IPLTF = 1 jIPLTL = 180 jNPLT = 1 

ITABL = l )IPROF = 1 

ABSIZ = 1.000E+00 50RSIZ = 1.000E+00 
NTZPRF = 1 JNCOLPF = 50 

ICHR = 2 iNCOLCO = AS jNROWCO = 20 

NO PATCHES YET USEO FOR THIS GROUP 



************************************************************ 
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GROUP 24. DUMPS FOR RESTARTS 
SAVE c T jAUTOPS = F jNOWIPE * F 
NSAVE =DF09 
LUSAVE = 9 

*** GRID-GEOMETRY INFORMATION *** 



X-COOROINATES 


OF THE CELL 


CENTRES 






1 . 550E-01 


S • 837E-01 


1 . 131E +00 


1 . 679E ♦00 


2 . 226E ♦OO 


3.500E+00 


5 . SOOE ♦00 


7 .500E+00 


9.200E+00 


1 . 060E +01 


Y -COORDINATES 


OF THE CELL 


CENTRES 






2.810E-01 


8.430E-01 


1.40SE+00 


1.843E+00 


2 . 15 7E +00 


2.595E+00 


3 . 157E+00 


3 . 719E+00 


4.375E+00 


S.12SE+00 


S.87SE+00 


6 . 62SE ♦00 








2-COORDINATES 


OF THE CELL 


CENTRES 






5 . OOOE-Ol 


1 . 850E +00 


3 . 700E+00 


5.550E+00 


7. 900E+00 


1 . 090E ♦01 


1.315E+01 


1 . 440E ♦01 


1.575E+01 


1. 710E+01 


1 . 810E +01 


1.945E^01 


2 . 080E +01 


2 . 180E ♦Ol 


2 . 258E+01 


2 . 315E+01 


2 . 372E>01 


2 . 417E ♦01 


2 . 452E ♦Ol 


2.620E+01 


2 . 920E+01 


3 . 220E+01 


3 .520E ♦Ol 


3 . 820E +01 


4. 120E+01 


4.420E+01 


4.670E+01 


4 . 807E+01 


4 . 832E+01 


4. 957E+01 


5. 032E +01 


5 . 099E+01 


S.1S6E+01 


S . 214E+01 


5 . 271E+01 


5.32S£^01 


5.375E+01 


5 .425S+01 


S . 475E ♦Ol 


5. 650E +01 



- INTEGRATION OF EQUATIONS BEGINS 



TIME STEP = 1 


SWEEP = 


2 








TOTAL 


RESIOUAL/l 


1.000E-06 ) 


FOR 


PI 


IS 


8.007E+05 


TOTAL 


RESIOUAL/l 


1.000E-06) 


FOR 


U1 


IS 


9. 092E*06 


TOTAL 


RESIOUAL/( 


1 . 000E-06 ) 


FOR 


VI 


IS 


1 . 282E+07 


TOTAL 


R£SIOUAL/( 


1.000E-06) 


FOR 


W1 


IS 


2 .376E+07 


TOTAL 


RESlOUAL/( 


1.000E-06) 


FOR 


KE 


IS 


1 . 210E +09 


TOTAL 


RESIDUAL/< 


1 . 000E-06 J 


FOR 


EP 


IS 


1.49cE+ll 


TOTAL 


RESIOUAL/( 


1.000E-06) 


FOR 


HI 


IS 


S . 642E ♦lO 


TIME STEP = 1 


SWEEP = 


3 








TOTAL 


RESlOUAL/( 


1 . 000E-06 ) 


FOR 


PI 


IS 


3.733E+05 


TOTAL 


RES10UAL/( 


1.000E-06) 


FOR 


U1 


IS 


1.050E+07 


TOTAL 


RESIDUAL/( 


1.000E-06) 


FOR 


VI 


IS 


1 . 474E ^07 


TOTAL 


RESIDUAL/1 


1.000E-06) 


FOR 


HI 


IS 


2.037E+07 


TOTAL 


RESIDUAL/! 


1.0C0E-06) 


FOR 


KE 


IS 


1.424E+09 



E3 



TOTAL RESIDUAL/I 
TOTAL RESI0UAL/( 
TIME STEP = 1 
TOTAL RESI0UAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/l 
TOTAL RESlOUAL/( 
TOTAL RESIDUAL/l 
TOTAL RESIOUAL/l 
TOTAL RESI0UAL/( 
TIME STEP = 1 
TOTAL RESI0UAL/( 
TOTAL RESI0UAL/( 
TOTAL RESIOUAL/t 
TOTAL RESI0UAL/( 
TOTAL RESIOUAL/l 
TOTAL RESI0UAL/( 
TOTAL RESIOUAL/l 
TIKE STEP = 1 
TOTAL RESIOUAL/l 
TOTAL RESIOUAL/l 
TOTAL RESI0UAL/( 
TOTAL RESI0UAL/( 
TOTAL RESIOUAL/l 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/l 
TIME STEP = 1 
TOTAL RESIDUAL/! 
TOTAL RESIOUAL/t 
TOTAL RESIOUAL/t 
TOTAL RESIOUAL/t 
TOTAL RESIOUAL/t 
TOTAL RESIOUAL/t 
TOTAL RESIOUAL/t 
TIME STEP = 1 
TOTAL RESIOUAL/l 
TOTAL RESIDUAL/t 
TOTAL RESIOUAL/t 
TOTAL RESIOUAL/t 
TOTAL RESIOUAL/t 



1 . 000E-06 ) 


FOR 


EP 


1, 000E-06 ) 


FOR 


HI 


SWEEP = 


174 




1. 000E-06) 


FOR 


PI 


1, 000E-06 ) 


FOR 


U1 


1. 000E-06 ) 


FOR 


VI 


1 . 000E-06 ) 


FOR 


HI 


1 . 000E-06 ) 


FOR 


KE 


1 . 000E-06 ) 


FOR 


EP 


1. 000E-06 ) 


FOR 


HI 


SHEEP = 


175 




1 , OOOE-06 ) 


FOR 


PI 


1. 000E-06 ) 


FOR 


U1 


1. OOOE-06 ) 


FOR 


VI 


1. 000E-06 ) 


FOR 


HI 


1. 000E-06 ) 


FOR 


KE 


1. OOOE-06 ) 


FOR 


EP 


1. OOOE-06 ) 


FOR 


HI 


SWEEP s 


176 




1 . 000E-06 ) 


FOR 


PI 


1. 000E-06) 


FOR 


U1 


1. 000E-06) 


FOR 


VI 


1. 000E-06) 


FOR 


HI 


1. 000E-06) 


FOR 


KE 


1. 000E-06 ) 


FOR 


EP 


1. 000E-06 ) 


FOR 


HI 


SWEEP = 


177 




1. 000E-06) 


FOR 


PI 


1. 000E-06) 


FOR 


U1 


1. 000E-06 ) 


FOR 


VI 


1. OOOE-06) 


FOR 


HI 


1. OOOE-06 ) 


FOR 


KE 


1. 000E-06 ) 


FOR 


EP 


1. 000E-06) 


FOR 


HI 


SHEEP = 


178 




1. 000E-06) 


FOR 


PI 


1. 000E-06) 


FOR 


U1 


1. OOOE-06 ) 


FOR 


VI 


1. 000E-06 ) 


FOR 


HI 


1. OOOE-06) 


FOR 


KE 



IS 2.341E+10 
IS 2 .620E+10 

IS 3 . 231E+04 
IS 6.021E+06 
IS 6 . 189E^06 
IS 8. 271E+06 
IS 2 . 926E *08 
IS 2 . 455E+10 
IS 2.619E+10 

IS 3 . 208E *04 
IS 5. 996E+06 
IS 6 . 183E *06 
IS 8. 329E*06 
IS 2 . 915E+08 
IS 2.393EtlO 
IS 2 . 628E+10 

IS 3.193E+04 
IS 6.071E+06 
IS 6 . 125E*06 
IS 8 . 259E*06 
IS 2.909E+08 
IS 2.4S5E+10 
IS 2 . 628E+10 

IS 3.235E+04 
IS 6.038E+06 
IS 6 . 128E+06 
IS 8.291E+06 
IS 2 . 934Et08 
IS 2.466E+10 
IS 2.622E+10 

IS 3.272E+04 
IS 6 . 095E+06 
IS 6 . 115E«-06 
IS 0 . 309E *06 
IS 3.010Et08 
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TOTAL RESIOUAL/< 1.000E-06) FOR EP IS 2.607E*10 
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3.115E 



FIELD VALUES OF U1 

IY= 12 1.596E+00 3.508E*00 4.603E+00 5.08SE+00 5.110E+00 ‘♦.405E*00 3.539E+00 2.288E+00 4.752E-08 
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FIELD VALUES OF VIST 

IY= 12 1.240E-06 1.173E-06 9.411E-07 6.256E-07 3.178E-07 8.973E-08 5.535E-08 1.576E-08 8.094E-11 9.000E-12 



IY= 11 7.568E-06 7.340E-06 6.560E-06 5.387E-06 3.952E-06 1.741E-06 9.548E-07 1.364E-07 2.682E-09 9.000E-12 
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FLOW FIELD AT ITHYD= 1, IZ= 31, ISWEEP= 180, ISTEP= 
FIELD VALUES OF PI 

IY= 12 -6.760E400 -1.177E+01 2.995E*00 2.855E+01 
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TOTAL RESIOUAL/I 1.000E-06) FOR PI IS 3.l85E*0<f 

TOTAL RESIOUAL/I 1.000E-06) FOR U1 IS 6.088E+06 



TOTAL RESIDUAL/* 1.000E-06) FOR VI IS 6 

TOTAL RESIDUAL/C 1.000E-06) FOR HI IS 8 

TOTAL RESIDUAL/* 1.000E-06) FOR KE IS 2 

TOTAL RESIDUAL/* 1.000E-06) FOR EP IS 2 

TOTAL RESIDUAL/* 1.000E-06) FOR HI IS 2 

WHOLE-FIELD RESIDUALS BEFORE SOLUTIONS 
WHOLE-FIELD SUM OF ABSIVOL.FLOW RESIDUALS) 
WHOLE-FIELD SUM OF ABS* RESIDUALS ) OF U1 
WHOLE-FIELD SUM OF ABS* RESIDUALS ) OF VI 
WHOLE-FIELD SUM OF ABS* RESIDUALS ) OF HI 
WHOLE-FIELD SUM OF ABS* RESIDUALS ) OF KE 
WHOLE-FIELD SUM OF ABS* RESIDUALS ) OF EP 
WHOLE-FIELD SUM OF ABS* RESIDUALS ) OF HI 
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NORMAL STOP REACHED IN PROGRAM 



DIRECTION | OF 



APPENDIX C 



AIRFLOW CHARACTERIZATION SHEET 



TEST PROCHAWi Testing of Navy Aircraft in Arched Structure Hush House 
LOCATION! Kelly A.F.B 
OaTCs _^JJovember_198^ — __ _ 

AIRCRAFT: F-4N 

ENQIHEi J-79-CE-8 Cel ‘- DEPRESSION (U. it.)! 1.7 

EVENT: 5 

CONDITIONS: Port Engine - A/B Pouer, Starboard Engine - Idle Pouer 
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fliRFi nw CHARACTERIZATION SHEET 



Thesis 

Iocat ion? RAM NavaT~Air Statio n. Jacksonville, Florida 
N/A 



N/A 



DATE; 

AIRCRAFT: , 

engines J-79-GE-8 

N/A 



CELL DEPRESSION (I*. »«■>.* N/A 



EVENT: 



CONDITIONS 



q-jnnlfi Engine with A/B @ centerline. 



AUOUENTER-INLET 
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APPENDIX D 



EXHAUST FLOW DATA SHEET 



yggl PROGRAM: Testing of Navy Aircraft Id Arched Structure Hush House 

LOCATION: ICelly A. F.B., San ftntopio, Texas 

DATE: 9 November 1982 

AIRCRAFT: F-4-N 

ENGINE: J-79-CE-8 

EVENT: •> 

CONDITIONS: Fort Engine - A/B, Starboard Engine - Idle Power 



AUGMENTED TUBE 




Tavg= L+5L f 

AUGMENTATION ratios 6.3 
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EXHAUST FLOW DATA SHEET 



TEST PROGRAM: Thesis 

location: Naval Air Station, 



DATE: 



AIRCRAFT: 



N/A 

n7a" 



ENGINE: J-79-GE-8 
EVENT: N/'k 



Jacksonv ille. Florida 



conditions: Single Engine with A/B @ centerline 



AUGMENTER TUBE 




FLOW (!&«./••*.)« 2 373 Vavg- 2 35.5 

augmentation ratio* 10. 6 Tavg= 401.6 
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APPENDIX E 

DATA TRANSFER BY TAPE 

The IBM 3033 main-frame presently does not have "PHOTON" 
graphics capabilities. In order to produce various visual 
plots that enhance the interpretation of tabulated data, the 
"DF09 " and "DF12" files must be properly formatted and then 
transferred to tape (for portability) allowing use of "PHOTON" 
on another computer which is "PHOTON" capable. 

The computer center will provide a non-labelled tape for 
data transfer. The operator, when queried by the programs 
below, will mount the tape allowing the "DF09" and "DF12" 
files to be transferred to the tape. 

To prepare the files for transfer to tape a copy is made 
to ensure that the original data file is safe. Program "FORM09 
JCL" properly formats and then writes the "DF09" file to a new 
file called "DF09C". Program "F0RM12 JCL" properly formats and 
then writes the "DF12" file to a new file called "DF12C". The 
data sets "DF09C" and "DF12C" which are produced using MVSHELP 
must be at least as large as the original "DF09" and "DF12" 
data sets and created before running "F0RM12 JCL" or "FORM09 
JCL" . 

With the two files now properly formatted the file "DSTP 
JCL" is used to transfer the designated file to tape. In the 
following example the "DF09C" file is transferred as the 
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second file to the tape called "NICOLA" . File "TSCAN JCL" then 
scans the tape to ensure the program was actually transferred 
to the tape. 

Two other files, "TAP JCL" and "DUMP JCL", can dump a 
file from tape to the virtual machine and return the first 
block of the tape file respectively. With the files on tape 
and checked the tape can be dismounted by the operator at the 
computer center and transported to a "PHOTON" capable system 
in the "formatted" form. 

The various files referenced are now listed. 



FILE: DSTP JCL A1 



/ /DSTP JOB ( 1541 ,9999 )> 'ERIC ' ,CLASS=E 
// EXEC PGM=IEBGENER 
//SYSPRINT DD SYSOUT=A 
//SYSIN OD DUMMY 

//SYSUT1 DD DISP=SHR ,DSN=MSS.S1541 . PHOENICS. DF09C 
//SYSUT2 DD UNIT=3400-6 , VOL =SER=NICOLA,DISP=( NEW, PASS ) , 

// LABEL = (2,NL,,0UT), 

// DCB = ( RECFM=FB,LRECL=80 ,BLKSIZE=800 >DEN=3 ,OPTCD=Q ) 
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FILE: TAP 



JCL 



A1 



//TAP JOB i 1541, 9999 ),*TAP JCL*,CLA$S=E 
//*MAIN CARDS=( 50 ) 

// EXEC PGM=IEBGENER 
//SYSPRINT DD SYSOUT=A 
//SYSIN DO DUMMY 

//SYSUT1 DD UNIT=3AOO-6,VOL=SER=NICOLA,DISP=( OLD >PASS ) , 
// LABEL=(2,NL,,IN), 

// DCB=(RECFM=FB,LRECL=80,BLKSIZE=800,OPTCD=Q) 

//SYSUT2 DD $YSOUT=B,DCB=BLKSIZE=80 
// 



FILE: DUMP JCL A1 



//DUMP JOB i 1541 >9999 ) , ' DUMP ' ,CLASS=E 
// EXEC TAPE ,VOLIN=NICOLA 
//SYSIN DD * 

DMPRECt 10,272,1 ),SKPFIL< 10 ,1 ) ,XECITTE( 999) 
/* 

// 



FILE: TSCAN JCL A1 



//GSCAN JOB ( 1541 ,9999 ) , 'NICOLA ' ,CLASS=E 
// EXEC TSCAN , VOLIN=NICOLA 
// 
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FILE: FORM09 JCL 



A1 



//FORM09 JOB (1591,9999), 'CHAM' >CLASS=B 
//*MAIN LINES = ( 20 ) 

// EXEC FORTVCLG 
//FORT.SYSIN DD * 

C PROGRAM FORM09 
C 

C Increase PARAMETER settings if more variables or a larger grid 
C is required. 

C 

PARAMETER! NPHIP=80 ) 

C MB NXYMAX must be >= NZ and NX*NY 
PARAMETER! NXYMAX=20000 ) 

CHARACTER** NAME ( NPHIP ), MESS! 10 ) 

CHARACTER*15 FILIN, FILOUT 

CHARACTER*1 ANS 

CHARACTER*10 ACCESS 

LOGICAL LBUF ( 10 ), STORE! NPHIP ), DEBUG 

DIMENSION RBUF! NXYMAX ),IBUF( 10) 

C 

100 FORMAT! IX, 71 10) 

101 FORMAT! 6(1 PE 13. 6)) 

102 FORMAT! IX, 79L1) 

103 FORMAT! IX, 19A9) 

C 

C Reset DEBUG=.TRUE. to get error messages.. 

DEBUG= .TRUE . 

C 

C Set Logical Units... 

LUIN=9 

LU0UT=10 

LUBUG=6 

C 

C Set File Names... 

C Name of unformatted restart file 
C WRITE ( 6 ,* ) 1 RESTART FILENAME? 1 

C READ! 5, ' ( A ) ' ) FILIN 
C Name of formatted restart file 
C WRITE! 6,*) ' FORMATTED FILENAME?' 

C READ! 5,' (A)' ) FILOUT 

C 

C Direct access or sequential? 

C WRITE! 6,*) ' DIRECT ACCESS (PHIDA) FILE? Y/N ' 

C READ! 5, '! A1 ) ' ) ANS 

C IF! ANS.EQ. 'Y ' .OR. ANS.EQ. 'Y' ) THEN 

C Version 1.9 default 
C ACCESS=* DIRECT' 

C LENREC=256 

C ELSE 

C Version 1.3 

ACCESS= 'SEQUENTIAL' 
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oo ooooooooooo 



FILE 



F0RM09 JCL 



A1 



LENREC=100 

ENDIF 

WRITE ( 6 ,* ) 1 DEFAULT LENREC? (100 VI. 3, 256 VI. 4) - Y/N' 

READ( 5, '(A)' ) ANS 

I F ( ANS . NE . * Y ' . AND . ANS . NE . ' Y ' ) THEN 
WRITE! 6,*) ’ LENREC? (WORDS NOT BYTES) 1 
READ! 5, ' ( I ) ' ) LENREC 
ENDIF 

NBYTES is the number of bytes per single precision word. 

NBYTES=4 

Record length in bytes. 

LREC=NBYTES# LENREC 

OPEN input file: 

OPEN! UNI T= 9, FORM= 'UNFORMATTED ' , STATUS 2 ‘OLD * , 

1 ACCESS=ACCESS,IOSTAT=IERR,ERR=90 ) 

IF! IERR.NE .0 ) GO TO 90 
IF! ACCESS. EQ. ’DIRECT' ) THEN 
READ! LUIN ,REC=1 ) ( MESS! I ) , 1 = 1 , 10 ) 

ELSE 

READ! LUIN) ( MESS! I ) ,1 =1 , 10 ) 

ENDIF 

IF! DEBUG) WRITE! LUBUG,# ) MESS 
IF! ACCESS. EQ. ’DIRECT’ ) THEN 
READ! LUIN, REC = 2) ( LBUF! I ) ,1 = 1,4 ) 

ELSE 

READ! LUIN ) ( LBUF ! I ) ,1 = 1 >4 ) 

ENDIF 

IF ! DEBUG ) WRITE t LUBUG,# ) ( LBUF! I ),I=1,4) 

IF! ACCESS. EQ. 'DIRECT* ) THEN 

READ! LUIN ,REC=3 ) NX, NY ,NZ,NPHI ,! IBUF! I ) ,1 = 1 ,6 ) , LENREC 
ELSE 

READ! LUIN) NX, NY ,N2,NPHI , ( IBUF! I ) ,1 = 1 ,6 ) , LENREC 
ENDIF 

NXNY=NX*NY 

I F ! NXNY . GT . NXYMAX ) GOTO 92 
I F ! NPHI . GT . NPHIP ) GO TO 93 
IF! DEBUG) WRITE! LUBUG,*) NX,NY,NZ 

OPEN output file: FORMATTED, SEQUENTIAL and of FIXED RECORD LENGTH. 
OPEN! UNIT=10 ,FORM= ’ FORMATTED ' ,STATUS= ’OLD ’ , 

1 ACCESS= ’ SEQUENTIAL ' , IOSTAT =IERR , ERR = 91 ) 

IF(IERR.NE.O) GO TO 91 
WRITE! LUOUT, 103 ) ! MESS! I ) ,1 = 1 , 10 ) 

WRITE! LUOUT, 102) ! LBUF! I ) ,1 = 1,4 ) 

WRITE! LUOUT, 100 ) NX, NY ,NZ,NPHI , ( IBUF! I ) ,1 = 1,6 ), LENREC 
IF! ACCESS. EQ. ’DIRECT' ) THEN 
READ! LUIN,REC=4 ) RINNER 
ELSE 

READ! LUIN) RINNER 
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FILE: FORM09 JCL 



A1 



ENDIF 

WRITE! LUOUT, 101) RINNER 

IF! ACCESS. EQ. 'DIRECT' ) THEN 

READ! LUIN, REC=5) (NAME! I ),I=1,NPHI ) 

ELSE 

READ! LUIN ) I NAME! I ) ,1=1 , NPHI ) 

ENDIF 

C Set C36... names for variables 51 to NPHI for version 1.4, to avoid 
C control characters etc I I Unless already stored) 

IF( ACCESS. EQ. 'DIRECT' .AND. NPHI .GT. 50 ) THEN 
READ! LUIN, REC = 10) I STORE! I ) ,1=1, NPHI ) 

DO 6 I =51, NPHI 
J=I-15 

IF! .NOT. STORE! I ) ) THEN 
IF! I. LT. 100) THEN 

WRITE! NAME (I ),'! "C" #12,* ' ")') J 
ELSE 

WRITE! NAME! I ),'( "C">13)') J 
ENDIF 
ENDIF 

6 CONTINUE 
ENDIF 

WRITE! LUOUT ,103 ) I NAME! I ), 1=1, NPHI ) 

IF! DEBUG) WRITE! LUBUG,*) ( NAME ( I ) ,1 = 1 ,NPHI ) 

IF(ACCESS.EQ. 'DIRECT' ) THEN 
READ! LUIN, REC = 6) ( RBUF! I ),I = 1,NX ) 

ELSE 

READ! LUIN ) t RBUF! I ) ,1=1 ,NX ) 

ENDIF 

WRITE! LUOUT, 101) I RBUF! I ) ,1 = 1, NX ) 

IF! ACCESS. EQ. 'DIRECT' ) THEN 
READ! LUIN, REC = 7) ( RBUF! I ) ,1 = 1, NY ) 

ELSE 

READ! LUIN) ! RBUF! I ) ,1 =1 ,NY ) 

ENDIF 

WRITE! LUOUT, 101) ! RBUF! I ) ,1 = 1 ,NY ) 

IF! ACCESS. EQ. 'DIRECT' ) THEN 
READ! LUIN, REC=8) ! RBUF! I ) ,1 =1,N2 ) 

ELSE 

READ! LUIN ) ! RBUF! I ),I=1,N2 ) 

ENDIF 

WRITE! LUOUT, 101) ! RBUF! I ) ,I = 1,N2 ) 

IF! ACCESS. EQ. 'DIRECT' ) THEN 
READ! LUIN, REC=9) ! RBUF! I ),I=1,N2) 

ELSE 

READ! LUIN ) ! RBUF! I ) ,I=1,N2 ) 

ENDIF 

WRITE! LUOUT, 101) I RBUF! I ),I=1,NZ) 

IFUCCESS.EQ. 'DIRECT' ) THEN 

READ! LUIN,REC=10 ) ( STORE ! I ) ,1 = 1 , NPHI ) 
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FILE: FORM09 JCL A1 



ELSE 

READ(LUIN) ( STORE ( I ) ,1 = 1 ,NPHI ) 

ENDIF 

WRITE< LUOUT >102 ) ( STORE ( I ) ,1 = 1 ,NPHI ) 

IF( DEBUG ) WRITE (LUBUG, 102) ( STORE ( I ) >1 = 1 >NPHI ) 

NVAR=0 

DO 5 I=1,NPHI 
IF(STORE(I)) NVAR=NVAR+1 
5 CONTINUE 
IREC=10 
DO 10 K=1,NZ 
DO 20 I=1,NVAR 
IF(ACCESS.EQ. ’DIRECT* ) THEN 
IREC=IREC+1 
NVREC=NXNY/LENREC 

IF(MOD(NXNY,LENREC).NE.O) NVREC=NVREC*1 

DO 30 IVREC=1,NVREC 

JF = ( IVREC-1 )*LENREC ♦ 1 

JL=IVREC*LENREC 

JL=MIN0( JL,NXNY) 

READ( LUIN,REC=IREC) ( RBUF( J ) , J=JF , JL ) 

30 cormNUE 
ELSE 

READ(LUIN) (RBUF( J),J = 1,NXNY) 

ENDIF 

WRITE (LUOUT, 101) ( RBUF ( J ) > J=1 ,NXNY ) 

20 CONTINUE 
10 CONTINUE 

IF( DEBUG) WRITE (LUBUG, 1900) 

1900 FORMAT ( IX, * NORMAL STOP IN PROGRAM') 

GO TO 999 

C Error trapping... 

90 IF( DEBUG) WRITE ( LUBUG, 1901 ) 

1901 FORMAT ( IX, * ERROR IN OPENING INPUT FILE') 

GO TO 999 

91 IF( DEBUG) WRITE ( LUBUG, 1911 ) 

1911 FORMAT( IX, 'ERROR IN OPENING OUTPUT FILE') 

GO TO 999 

92 I F ( NXNY . GT . NXYMAX ) WRITE ( LUBUG, 1921 ) 

I F ( N2 . GT . NXYMAX ) WRITE ( LUBUG , 1923 ) 

1921 FORMAT( IX,' SETTING OF NX*NY TOO LARGE) RESET IN FORM09 AND * 
1 ,/, IX, 'RE-CREATE THE LOAD MODULE*) 

1923 FORMAT( IX, 'SETTING OF NZ TOO LARGE) RESET IN FORM09 AND * 
1,/, IX, 'RE-CREATE THE LOAD MODULE') 

GO TO 999 

93 WRITE ( LUBUG ,193) 

193 FORMAT( IX, 'SETTING OF NPHI TOO LARGE) RESET IN FORM09 AND ' 
1,/, IX, 'RE-CREATE THE LOAD MODULE') 

999 CONTINUE 
C CLOSE files ... 
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CLOSE ( LUIN > STATUS =* KEEP ' ) 

CLOSE(LUOUT,STATUS='KEEP* ) 

STOP 

END 

//GO. FT09F001 DD DSN=MSS.S1541 .PHOENICS.DF09,DISP=(OLD ) 
//GO. FT10F001 DD DSN=MSS. S1541 . PHOENICS. DF09C ,DISP = ( OLD ) 
/* 

// 



137 



FILE: F0RM12 JCL 



AI 



//F0RM12 JOB (1541, 9999 ),*CHAM',CLASS=B 
//*MAIN LINES=( 20 ) 

// EXEC FORTVCLG 
//FORT.SYSIN DD * 

C File name: FORM12 Iasi amended 15 October 1986 
C PROGRAM FORM12 

C Program to rewrite a DF12 file into a formatted READCO file. 

C (JPE & JFD ) 

C 

C Reset NXY1P if larger grid required. . 

PARAMETER( NXY IP =10000 ) 

DIMENSION AREC(NXYIP) 

C Reset fixed record length of DF12 file if required.... 

CHARACTERS FILIN, FILOUT 
LOGICAL DEBUG 
INTEGER SECTOR 
DATA LENREC/100/ 

C 

100 FORMAT ( 315 ) 

101 F0RMAT(5( 1PE13.6)) 

C 

C SECTOR represents the largest record length that can be OPENed. 
SECT0R=65536 

C NBYTES is the number of bytes per single precision word. 

NBYTES=4 

LREC=MIN0( NBYTES*LENREC,SECTOR-NBYTES ) 

IBLK=LREC+NBYTES 

C 

C Set Logical Units... 

LUIN=12 
LUOUT = 13 
LUBUG=6 

C Reset DEBUG if error messages required.. 

DEBUG= . TRUE . 

C Set file names 
C FI LIN= 1 DF 12 

C FILOUT = 1 DF 12 FORM 1 

C 

C OPEN input file: DIRECT ACCESS, UNFORMATTED and of FIXED 
C RECORD LENGTH. 

C 

OPENt UNIT=12,F0RM= 'UNFORMATTED ' ,STATUS= 'OLD ' , 

1 ACCESS= ' DIRECT ' ,RECL=LREC , 

1 IOSTAT =IERR,ERR=90 J 

IF(IERR.NE.O) GO TO 90 
C 

READ( LUIN ,REC = 1 ) NXP1 ,NYP1 ,N2P1 
NXY1=NXP1*NYP1 

IF ( NXY1 .GT.NXY1P ) GO TO 92 
NREC=NXY1/LENREC 
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ISH0RT=NXY1-LENREC*NREC 

IF! ISHORT . GT . 0 ) NREC=NREC+1 
IF! ISHORT . EQ. 0 ) ISHORT=LENREC 

IF! DEBUG) WRITE! LUBUG,*) NXP1 >NYP1,NZP1>NXY1,NREC , ISHORT 

OPEN output file: FORMATTED, SEQUENTIAL and of FIXED RECORD LENGTH. 

OPEN! UNIT=13 , FORh= ' FORMATTED ' ,STATUS= ' OLD ' , 

1 ACCESS= 'SEQUENTIAL ' ,IOSTAT=IERR,ERR = 91 ) 

IF! IERR.NE. 0 ) GO TO 91 
WRITE! LUOUT ,100 )NXP1 ,NYP1 >NZP1 
NR = 1 

DO 5 K=1,N2P1 
DO 6 IC=1,3 
IMAX=LENREC 

DO 7 JR=1,NREC 

IF! JR.EQ.NREC ) IMAX=ISHORT 
NR=NR+1 

NRA=( JR-1 )*LENREC 

READ! LUIN,REC=NR) I AREC! NRA+I ) ,1=1 ,IMAX ) 

7 CONTINUE 
C Write data to output file 

WRITE! LUOUT, 101 )( AREC! II ) , 11 = 1 ,NXY1 ) 

6 CONTINUE 
5 CONTINUE 

IF! DEBUG) WRITE! LUBUG, 1900) 

1900 FORMAT! IX, 'NORMAL STOP IN PROGRAM') 

GO TO 999 

C Error trapping. . . 

90 IF! DEBUG) WRITE! LUBUG , 1901 ) 

1901 FORMAT! IX, 'ERROR IN OPENING INPUT FILE') 

GO TO 999 

91 IF! DEBUG) WRITE! LUBUG ,1911 ) 

1911 FORMAT! IX, 'ERROR IN OPENING OUTPUT FILE') 

GO TO 999 

92 IF! DEBUG) WRITE! LUBUG ,1921 ) 

1921 FORMAT! IX, 'GRID TOO LARGE > RESET NXY1P AND RE-CREATE', 

1 /, IX, ' LOAD MODULE' ) 

999 CONTINUE 
C CLOSE files . . . 

CLOSE! LUIN,STATUS= 'KEEP' ) 

CLOSE! LUOUT, STATUS= 'KEEP ' ) 

STOP 

END 

/ /GO . FT12F001 DD DSN=MSS . S1541 . PHOENICS . DF12 ,DISP=( OLD ) 

//GO . FT13F001 DD DSN=MSS.S1541 .PHOENICS. DF12C,DISP=< OLD ) 

/* 

// 
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